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We review results from lattice QCD calculations on the thermodynamics of strong- 
interaction matter with emphasis on input these calculations can provide to the 
exploration of the phase diagram and properties of hot and dense matter created 
in heavy ion experiments. This review is organized as follows: 

Introduction 

QCD thermodynamics on the lattice 
QCD phase diagram at high temperature 
Bulk thermodynamics 
Fluctuations of conserved charges 
Transport properties 
Open heavy flavors and heavy quarkonia 
QCD in external magnetic helds 
Summary 


1. Introduction 

It has long been recognized that under extreme conditions of high temperature or 
densities matter interacting through the strong force - strong-interaction matter 
- cannot exist in the form of nuclear matter formed from hadrons. The copious 
production of new resonance^ and the intrinsic size of the nucleons put a natural 


1 










2 


H.-T. Ding, F. Karsch and S. Mukherjee 


limit on the range of validity of hadron physics EEl In the context of thermodynam¬ 
ics of strong-interaction matter this found its expression in the formulation of the 
hadron resonance gas modellEI which led to the postulate of a limiting temperature 
- Hagedorn temperature - for the thermodynamics of ordinary nuclear matter. 

With the formulation of Quantum Chromodynamics (QCDjUH as the funda¬ 
mental theory of strong interactions it has soon been realized that matter formed 
from strongly interacting particles may be converted into a new form of matteiEISl 
-the quark-gluon plasmgP - in which the dominant degrees of freedom are the con¬ 
stituents of hadrons, i.e. quarks and gluons. At low temperatures and densities 
they are confined in hadrons. However, at high temperature and densities they can 
move “freely” over macroscopic distances in bulk strong-interaction matter. 

It has been speculated that the transition from confined to deconfined matter, 
like in condensed matter systems, goes alon^I^l with the restoration of global sym¬ 
metries of the underlying microscopic theory describing the interaction among con¬ 
stituents of strong-interaction matter, i.e. QCD. The analysis of these symmetries 
led to the phenomenology of a complex phase structure of strong-interaction mat¬ 
ter, leading to homogeneous and inhomogeneous as well as color-superconducting 
phases. This bears many analogies to phase structures known in condensed matter 



Fig. 1. QCD phase diagram in the temperature (T) and baryon chemical potential {i^b) plane. 
At vanishing iib lattice QCD calculations show that the transition is not a phase transition but a 
continuous crossover, reflecting pseudo-critical behavior in the vicinity of the true chiral transition. 
At fiB > 0 a second order critical end point (CEP) may exist, which would be followed by a line 
of first order transitions at larger values of fiB- Constraints on the location of a CEP will come 
from lattice QCD calculations. 
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systems EEH Only little of this complex structure is verified through first principle 
lattice QCD calculations. In Fig.[^we show a sketch of the QCD phase diagram that 
only highlights those aspects that we believe can be addressed at present through 
lattice QCD calculations. We will discuss the phase diagram in detail in Section 

Phase transitions as well as properties of matter close to such transitions arise 
from complex long range interactions in multi-particle systems on multiple length 
scales and at the same time often develop simple universal behavior. Their anal¬ 
ysis requires non-perturbative techniques. With the formulation of lattice QCD, 
i.e. the introduction of a discrete space-time lattice as regulator for quantum field 
theories like QCEP^a framework for new non-perturbative approaches to the study 
of strong-interaction matter had been established. It became possible to perform 
first-principle numerical calculations - lattice QCD simulationJI^I _ py which a sys¬ 
tematic study of the phase structure and basic bulk properties of strong-interaction 
mattei^i^K^ could be performed for the first time. We will give a brief introduction 
to lattice QCD in Section Today these calculations are highly advanced and con¬ 
tinuum extrapolated results for the QCD transition temperature and the equation 
of state exist at vanishing net baryon number density. We will review these results 
and report on the status of the lattice QCD program that aims at an extension of 
these calculations to non-zero net baryon number density or equivalently non-zero 
baryon chemical potential in Section 

Properties of strong-interaction matter are studied intensively in large exper¬ 
imental programs at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
National Laboratory (BNL), USA, and the Large Hadron Collider (LHC) as well 
as the Super Proton Synchrotron (SPS) at the European Research Center (CERN). 
Experiments at these accelerators are devoted to the exploration of the properties of 
hot and dense matter created in collisions of ultra-relativistic heavy ions. Through 
the variation of beam energies at RHIC, i.e. the beam energy scan (BES), it also 
became possible to explore systematically the phase structure of strong-interaction 
matter at non-zero net baryon number density. This allows to search systematically 
for the chiral critical point or critical end point (CEP) in the QCD phase diagram. 
The CEP is a postulated second order phase transition point which is expected to 
mark the endpoint of a line of first order phase transitions that separates the low 
temperature, low density hadronic phase from a low temperature, large baryon num¬ 
ber quark-gluon plasma phase (also the existence of this first order line at present 
is not confirmed through lattice QCD calculations). The existence of a CEP in 
the QCD phase diagram is imprinted in the properties of fluctuations of conserved 
charges, e.g. net baryon number, electric charge or strangeness fluctuations. These 
conserved charge fluctuations can be calculated in lattice QCD at vanishing baryon 
chemical potential and the calculations can be extended to non-zero values of the 
chemical potential using Taylor expansions of thermodynamic observables. They 
can directly be compared to experimental measurements of conserved charge fluc¬ 
tuations performed at the LHC and the BES at RHIC and may provide insight into 
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the existence and location of a CEP in the QCD phase diagram. We will report in 
Section!^ on the status of lattice QCD calculations of conserved charge fluctuations 
and their usage to analyze freeze-out conditions in heavy ion experiments. 

Heavy ion experiments at RHIC provided striking evidence for the highly cor¬ 
related, non-perturbative structure of strong-interaction matter in the high tem¬ 
perature phase. In particular, the observation of strong elliptic flow and very ef¬ 
ficient energy loss of high momentum quarks traversing through the medium (jet 
quenching) paved the way for a picture of a strongly interacting medium that has 
been called an ‘almost perfect liquid’. These measurements suggest that strong- 
interaction matter at temperatures close to but above the transition temperature 
has quite unique transport properties, (i) a small shear viscosity to entropy ratio 
that may be close to the conformal limit value r]/s = l/47r calculated in conformal 
field theories (AdS/CFT limit), (ii) a large bulk viscosity that may diverge in the 
massless QCD limit at the chiral phase transition temperature, (iii) small electrical 
conductivity and diffusion coefficients. These transport properties can systemati¬ 
cally be analyzed in lattice QCD, although this requires the reconstruction of con¬ 
tinuous spectral functions from a finite set of observables, which is in general an 
ill-posed problem. The calculation of transport properties in lattice QCD will pro¬ 
vide crucial input to the modeling of the hydrodynamic expansion of hot and dense 
matter formed in heavy ion collisions. Spectral functions also carry information on 
in-medium modifications of hadrons. The dissociation of heavy quark bound states 
has long been advocated as one of the striking signatures for the formation of a 
deconfined color screened mediurrP^ while the melting of light quark bound states 
is intimately related with the restoration of chiral symmetry.^^^ We will report on 
the status of calculations of transport coefficients and in-medium hadron properties 
in Sections [6] and O 

Finally, in Sectionwe will briefly review a relatively new but very exciting field 
of QCD thermodynamics that addresses the behavior of strong-interaction matter in 
strong external magnetic fields. Strong, static magnetic fields will certainly influence 
the thermodynamics and lead to modifications of the transition temperature as well 
as the equation of state. Strong fields are generated in the early phase of heavy ion 
collision£2l2Il but weaken rapidly during the subsequent expansion phase. While 
lattice QCD calculations will allow to understand basic effects that occur in matter 
exposed to strong fields it is at present not clear to what extent equilibrium lattice 
QCD calculations can contribute to a quantitative analysis of experimental data. 


2. QCD thermodynamics on the lattice 

2.1. Path integral formulation of QCD thermodynamics 

The equilibrium thermodynamics of elementary particles interacting only through 
the strong force is controlled by the QCD partition function which can be expressed 
in terms of a Euclidean path integral. The grand canonical partition function. 
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Z{T^V^jl)^ is given as an integral over the fundamental quark (t/S, and gluon 
(Ajy) fields. In addition to its dependence on volume (V), temperature (T) and a 
set of Nf chemical potentials, jl = RcC •••) fhe partition function implicitly 
depends on the masses, m = {mu^md^rris^ of the Nf different quark flavors. In 
this review we will most of the time discuss the thermodynamics of QCD with 2 
light quarks {u,d) which are assumed to have degenerate masses, rUu = rrid, and a 
heavier strange quark ( 5 ), with mass m^. This often is referred to as (2+l)-flavor 
QCD, or = 2 + 1. 

In Euclidean space-time, which is obtained from the Minkowski formulation by 
substituting t -A- —it with r G M, the QCD Lagrangian is given by 

^QCD ^gluon ^f ermion 

=^ , ( 1 ) 

f=u,d-,S‘“ 

where Greek letters are spinor indices, a = 1,.., — 1 is the color index, Nc is the 

number of colors {Nc = 3 for QCD) and is the mass of quarks with flavor /. 
The covariant derivative Ip^ and the field strength tensor are given by 

= i^N = > ( 2 ) 

F;, = d,Al - - gC^A^Al . (3) 

Here are the gauge fields, {'^J) are the quark (anti-quark) fields, are 
the generators of SU(A^c )5 fabc are the corresponding structure constants, 7 ^ are 
the Euclidean Dirac matrices obeying { 7 ^, 7 ^^} = and g is the bare coupling 

constant. 

In the Euclidean path-integral formalism the partition function of QCD is then 
given by 

Z{T,V,p)= fl[VA^ n (4) 

/X f = U,cLs... 

with the Euclidean action 

l/T 

SE{T,V,il) = - JdxoJdAC^ifl). (5) 

0 y 

Here we have suppressed the dependence of the Euclidean Lagrangian and action on 
the fields but have stressed explicitly their dependence on the various 

quark chemical potentials that couple to the conserved quark number currents 

H{ji) = c^cD + ■ 

f=u,d-,S“ 


( 6 ) 
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The thermal expectation value of physical observables O can be obtained 
through 

= Z(T,V,/2) /• (7) 

M / 


Basic thermodynamic quantities like the pressure (P), energy density (e) or 
net-quark number density Uf can be obtained from the logarithm of the partition 
function using standard thermodynamic relations, 

T = ^lnZ(TV,M), (8) 

e _ 1 dlnZ{T,V,il) 

F ~ ~vF dl/T 

Uf _ 1 d\nZ{T,V,fl) 

F ~ vF ’ 

where we introduced the chemical potentials in units of temperature, fij = i^f/T. 
In fact, the QCD partition function depends on chemical potentials only through 
these dimensionless combinations, which are the logarithms of the fugacities, Zf = 
exp(p//T). 

As will become clear in the next subsection, the numerical analysis of thermody¬ 
namic observables at non-vanishing chemical potential, /2 7 ^ 0 , is difficult in lattice 
regularized QCD. A viable approach that circumvents the so-called sign-problem in 
lattice QCD at 7 ^ 0, suitable for moderate values of the chemical potentials, is to 
consider Taylor expansions for thermodynamic observables. The starting point for 
such an analysis is the Taylor expansion of the pressure. In ( 2 -M)-flavor QCD it 
reads 



where the expansion coefficients are dimensionless, generalized susceptibilities that 
can be evaluated at /i = 0 , 



xtFT) 


Qi+j+k p 


( 12 ) 


We will present results from Taylor expansions of QCD thermodynamics in several 
chapters of this review. Other thermodynamic observables and expressions more 
suitable for the calculation of P/T^ and e/T^ in lattice regularized QCD will be 
introduced in later sections where they appear. 


2.1.1. The high temperature, ideal gas limit 

At very high temperatures, all quark masses are small on the scale of the temper¬ 
ature and massless QCD becomes a good approximation. Furthermore, because 
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of asymptotic freedom also interactions among quarks and gluons become small at 
high temperature. QCD thermodynamics thus approaches that of a non-interacting, 
massless quark-gluon gas. In the temperature range accessible to heavy ion colli¬ 
sions the massless limit of 3-flavor QCD is most relevant {Nf =3). In this limit 
the pressure is given by 


(T) =^+ y + 

ideal 45 2 It; inAT), 


(13) 


where the first term gives the contribution of gluons and the sum yields the contri¬ 
bution of quarks with different flavor degrees of freedom. 

The thermodynamics of QCD at high temperatures can be systematically an¬ 
alyzed in perturbation theory. This, however, becomes complicated beyond 0{g‘^) 
due to the appearance of non-perturbative length scales of 0{gT) and 0{g‘^T) re¬ 
flecting electric and magnetic screening lengths,!^ which require the resummation 
of diagrams leading to the so-called hard thermal loop perturbation theorypHH 
or the explicit integration over hard scales leading to the dimensional reduction 
scheme mfeU 


2.1.2. The low temperature hadron resonanee gas approximation 

At low temperature quarks and gluons are confined in colorless hadrons, i.e. baryons 
and mesons are the relevant degrees of freedom. In fact, it turns out that even 
at temperatures close to the transition region from hadronic matter to the high 
temperature quark-gluon plasma phase a non-interacting gas constructed from all 
experimentally known resonances does provide quite a good description of QCD 
thermodynamics. This had been anticipated by R. Hagedorn when he formulated 
the hadron resonance gas (HRGPI model prior to QCD. 

In an HRG model the pressure and other thermodynamic observables are easily 
obtained from the logarithm of the partition function. 

In Zhrg{T, C m) = ^ In (T, V,il)+ ^ In (T, R, fl) , (14) 

iG mesons iG baryons 

where fl = {rb, Mq, Rs) is the set of baryon number, electric charge and strangeness 
chemical potentials, respectively. The partition functions for mesons (M) or baryons 
(B) are given by 

V 

lnZ^/^{T, V, p) = dk e ln(l T , (15) 

with energies + m? and fugacities 

Siiis)/T) . (16) 

Of course, = 0 for all mesons and Bi = Ml for baryons. The set of hadron 
chemical potentials, related to conserved quantum numbers, and the set of quark 
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flavor chemical potentials are easily related to each other, 


= 3MB + 3^0 




Ms 


Mb 


7Mb 


Mq 


:Mq - M5 


(17) 


With this it is straightforward to rewrite the Taylor series given in Eq. in terms 
of quark chemical potentials, also in terms of baryon number, electric charge and 
strangeness chemical potentials. 



where the expansion coefficient^ can again be evaluated at /I = 0, 


(18) 


(19) 


2.2. Lattice Regularization and Continuum Limit 

Lattice gauge theory suggested by K. G. Wilson in 1974^ is based on the Euclidean 
path integral formalism and provides a particular regularization scheme for QCD 
by introducing a finite lattice spacing as cut-off. In this way space-time is dis¬ 
cretized and the path integral becomes a finite, yet high dimensional integral over 
the gauge and fermion field variables. The discretization of space and time intro¬ 
duces a systematic cut-off dependence in all observables, which vanishes when the 
lattice spacing is taken to zero, i.e. in the continuum limit. Many theoretical and 
technical details can be found in excellent textbooks such as Ref. [27^301 Here we 
will only give a brief introduction that may help to make this review self-contained. 

By introducing a hyper-cubic lattice of size x with a small but finite lattice 
spacing a the calculation of observables given by Eq. reduces to the evaluation of 
finite, but high dimensional integrals. Rather than calculating the partition function 
directly one usually is concerned with the calculation of expectation values, Eq. 
This can be done by exploiting well-known Monte Carlo simulation techniques.^ 
In order to preserve gauge invariance of the discretized action one introduces the 
gauge degrees of freedom as variables living on links between neighboring sites of the 
lattice, Ux,ji = ex.p{iagA^{x)), while the fermionic degrees of freedom, and 
are defined on the sites of the lattice. The latter are anti-commuting Grassmann 
variables and the former are Nc x Nc matrices that are elements of the SU{Nc) 
color group. 


^We will in general suppress super- and subscripts of the generalized susceptibilities, if a subscript 

BQS _ BS BQS _ B 

IS zero, i.e. X102 = Xf2 or X200 = X2 ■ 
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The volume V and the temperature T are related to the spatial and temporal 
extents of the lattice, respectively, 

Here N^- and N^- are the number of sites in spatial and temporal directions. A 
consequence of introducing a discrete space-time lattice is that aside from Lorentz 
symmetry also some of the symmetries of QCD get explicitly broken by the finite 
lattice cut-off. They will be recovered in the continuum limit, a ^ 0 at fixed V and 
T. Asymptotically, for small values of the bare gauge coupling g this is controlled 
by the QCD /d-function 


aAL = 



bil2bl 

^-l/2bog^ ^ 


( 21 ) 


with bo = - iNf) /167r2 and bi = (f iV^ - (^fN, + TV/) 

Keeping physical observables constant when approaching the continuum limit 
also requires a proper tuning of the bare quark masses. In the discretized version 
of QCD (lattice QCD) all quark masses are naturally expressed in units of the 
lattice spacing, ifif = rrifa. Taking the continuum limit thus requires to take the 
limit, g‘^ ^ 0 and m/ ^ 0 while keeping some physical observables, i.e. one per 
flavor degree of freedom, constant. This defines “lines of constant physics” (LCP) 
in the space of quark masses and the gauge coupling. In the case of (2+l)-flavor 
QCD with degenerate light quark masses, rriu = rrid^ one thus can determine the 
bare quark mass parameters mi = ihu and the strange quark mass ttIs using two 
physical observables. An often used approach is to fix the strange quark mass using 
the (fictitious) Rss meson mass, which only contains strange quarks and is quite 
insensitive to light quark mass values. This mass is matched to the lowest order 
chiral perturbation theory estimate — m^. The light quark mass is 

then fixed using a constant ratio rhs/rhi^ the value corresponding to the physical 
pion mass value is rhs/rhi = 27.5.^^ 

The lattice regularized partition function may be written as 


Z{T,V,fl) 


I H dip^jdip^j e , 

X,fi X,f 


I 


Y[detMfinf) e , 

f 


( 22 ) 


where Sg and S / are the discretized versions of the gluonic and fermonic part of the 
action Se- As the action is bilinear in the quark fields these can be integrated out, 
which gives rise to the determinant of fermion matrix Mf. Quite often exploratory 
lattice QCD studies are performed in the so-called “quenched approximation”. This 
refers to the calculation of observables that contain fermionic degrees of freedom 
as valence quarks. However, back-reaction of fermions on the gauge fields (virtual 
quark loops) are not included in the observables. This is achieved by fixing det Mf = 
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const, in Eq.[^ i.e. the determinant does not contribute in the generation of gauge 
field configurations nor does it contribute in the calculation of expectation values. 

The fermion matrix depends on the quark chemical potential /i/ which is in¬ 
troduced on the time-like links of the lattice.!^ Although other formulations are 
possibl^^a common approach is to replace the link variables, ^ exp {/gfa)U^ q 
and ^ exp {—/gfa)U^ The fermion determinant is positive definite only for 
fif = 0. For non-zero values of the chemical potential detM/(/i/) is a complex 
function. Still the partition function is real, i.e. integration over the gauge fields 
eliminates the imaginary contributions. However, the real part of the fermion de¬ 
terminant still changes sign, which prohibits the application of conventional Monte 
Carlo simulation techniques. This is known as the ’sign-problem’ in finite density 
QCD. There are various attempts to circumvent this sign-problem in numerical 
simulations performed directly at non-vanishing ft. At present complex Langevin 
simulation techniqueand the integration over Lefschetz thimble^^^K^are most 
actively being explored. In this review we will focus on the Taylor expansion ap- 
proaclEEll for the evaluation of thermodynamic observables at non-zero chemical 
potential. This only requires numerical simulations at vanishing baryon chemical 
potential. However, it is naturally limited in applicability through the radius of 
convergence of the Taylor series which remains to be determined for QCD. 

2.3. Fermion discretization schemes 

When discretizing the Euclidean action one has quite some freedom. Discretization 
schemes that differ by higher order corrections in the cut-off are equivalent to the 
extent that they yield the same physical answers in the continuum limit. The ba¬ 
sic construction principle of a lattice discretized version of QCD is that one wants 
to start with a theory that is manifestly gauge invariant. In addition one tries 
to choose discretization schemes that preserve as many of the symmetries of QCD 
as possible already at non-zero values of the lattice spacing. The discretization of 
the gluonic part of the action, 5'^, is rather straightforward and well-established im¬ 
provement schemes (Symanzik improvemenfElHll) are known that allow to eliminate 
systematically higher order cut-off errors. 

Some difficulties already arise in the naive discretization of the fermionic part of 
the Euclidean action, Sf. When one simply approximates the covariant derivative 
in the fermion action by nearest-neighbor differences of the quark fields additional, 
unphysical degrees of freedom arise even in the continuum limit owing to the pe¬ 
riodicity of the fermion dispersion relation on the lattice. This is known as the 
fermion doubling problem which is closely related to the explicit breaking of chiral 
symmetry in most commonly used fermion discretization schemes 

There are several discretization schemes for the fermionic part of the action such 
as Wilson fermions,^ staggered fermionsDomain Wall fermionas well as 
overlap fermions EEIl They differ to the extent they preserve chiral symmetry 
and/or eliminate the problem of fermion doublers. Wilson fermions avoid the don- 
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bling problem by adding a dimension five operator to the naively discretized fermion 
action!^ A variant of this is the so-called clover-improved Wilson fermion action 
where lattice cutoff effects are reduced from 0{a) to Most commonly 

used in QCD thermodynamics calculations is the staggered fermion discretization 
scheme in which the Dirac 4-spinor is spread over several lattice sitesThis reduces 
the doubler problem but does not eliminate it completely. However, it preserves a 
U (l)even ^U{l)odd remnant of chiral symmetry, which allows independent rotations 
of the quark fields on even and odd sites of the lattice, respectively. This is of great 
advantage in the study of chiral symmetry breaking at non-zero temperature as it 
allows to introduce an order parameter, the chiral condensate, that is sensitive to the 
spontaneous breaking of this continuous symmetry. The reduction of the doubler 
problem in the staggered fermion approach leads to the introduction of additional 
heavier states for each fermion flavor. In the continuum limit these so-called taste 
symmetry partners will become degenerate with the Goldstone mode that corre¬ 
sponds to the broken U{l)even x U{l)odd symmetry. The Goldstone boson thus has 
another 15 heavier partners (tastes) which become degenerate with the Goldstone 
boson only in the continuum limit. To reduce the taste symmetry breaking, i.e. 
reduce the mass difference between physical states and their heavier taste partners 
two improved staggered actions, the Highly Improved Staggered Quarks (HISQ)I^ 
and the stout smeared action,!^ are commonly used in lattice studies of QCD ther¬ 
modynamics. In Fig. 1^ (left) we show some results for the cut-off dependence of 
the root-mean-square mass values of the Goldstone particle and its 15 taste 

partners calculated in different staggered fermion discretization schemes. As seen in 
the left panel of Fig. the 4stout and HISQ actions have the smallest The 

right hand panel of this figure shows the influence of cut-off effects on the calculation 
of the fermion contribution to the pressure in the high temperature ideal gas limit. 
These cut-off effects are independent of the taste improvement schemes and solely 
arise from the strategy used to discretize the covariant derivatives appearing in the 
fermion action. The stout-smeared actions utilize the naive l-link discretization 
scheme that leads to 0{o?) errors in the ideal gas limit. In the HISQ, asqtad and 
p4^ actions 3-hnk terms are added to the action that eliminate the leading order 
cut-off effects in the ideal gas limit. In these cases discretization errors only start 
at 0{a^). 

Domain Wall fermions (DWF) and overlap fermions allow for an exact repre¬ 
sentation of chiral symmetry even at non-zero lattice spacing. In the Domain Wall 
discretization scheme the physical chiral Dirac fermions are constructed on two 4- 
dimensional hyper-surfaces at the edges of a 5-dimensional lattice with extent Lg in 
that fifth direction. These hyper-surfaces represent the usual 4-dimensional space- 
time and on each of the surfaces fermions with one given chirality exist. There 
persists a small explicit breaking of chiral symmetry which vanishes exponentially 
in the limit Lg ^ oc. In DWF calculations this gives rise to a small additive renor¬ 
malization of the quark masses, the so-called residual mas^m^es- Gontrolling and 
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Fig. 2. (Left) Root mean squared pion mass as a function of lattice spacing a. RMS 

Mtt is defined as the rooted sum of the mass squared of 16 pseudo scalar states divided by 4. The 
lattice cutoff effects are smaller with smaller (Right) Ratio of quark contribution to the 

pressure, obtained from lattice QCD calculations in the infinite temperature, ideal gas limit, to the 
pressure of an ideal quark gas in the continuum (Pideal) cis a function of Nr- P/Pideal approaches 
unity in the continuum limit. 


reducing these residual mass effects is one of the important improvement steps in 
calculations with DWF.^^ Overlap fermions on the other hand preserve exact chiral 
symmetry on a 4-dimensional lattice by obeying the Ginsparg-Wilson relation at 
non-zero lattice spacing.l^ 

Numerical calculations with Domain Wall as well as overlap fermions are quite 
time consuming. However, with increasing speed of super-computers calculations 
with physical light quark masses become feasible and chiral fermions have been used 
recently also for QCD thermodynamics studies.l^^M^ In particular in the analysis of 
subtle aspects of the QCD transition related to the temperature dependence of the 
axial anomaly, calculations with fermions obeying exact chiral symmetry already at 
non-zero values of the cut-off are mandatory. 

3. QCD phase diagram at high temperature 

Our thinking about the phase structure of strong-interaction matter centers around 
two very basic concepts in strong interaction physics - confinement and chiral sym¬ 
metry breaking. The former expresses the fact that only colorless states, baryons 
and mesons, can exist in the vacuum and are observed experimentally. This gave 
rise to the concept of a linearly rising, confining potential exhibited between quarks 
and anti-quarks, 

Vqq{r) = + err , (23) 

with CF being the string tension and a{r) the running coupling of QCD. 

Chiral symmetry breaking, on the other hand is a mandatory feature of strong 
interactions needed to explain the appearance of a light, almost massless, particle 
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in the hadron spectrum - the pions. In the limit of vanishing quark masses the 
QCD Lagrangian has a build-in chiral symmetry. It is invariant under independent 
global rotations of quarks with left and right handed chirality in flavor space as well 
as chiral rotations of single flavor quark spinors. This gives rise to the x 

U{Nf)ji chiral symmetry which is equivalent to SU{Nf)L x SU{Nf)ji x U{1)a x 
I/ y(l). The latter (1) reflects the conservation of baryon number, and the axial 
U{1)a symmetry, although an exact symmetry of the classical QCD Lagrangian, is 
explicitly broken by quantum fluctuations. The global SU{Nf)L x SU{Nf)ji flavor 
symmetry, on the other hand, is broken spontaneously in massless QCD, giving rise 
to massless Goldstone modes and a non-vanishing chiral condensate. 

The masses of the light pseudo-scalar (Goldstone) pions are then understood as 
arising from the small non-zero values of the light up and down quark masses. This 
is reflected in the Gell-Mann-Oakes-Renner relation between the pion mass {m^), 
the quark masses {rriu^md), the pion decay constant and the non-vanishing 
chiral condensate = {ip'ip)u + ('0'0)d fhat arises from the spontaneous breaking 
of chiral symmetry, 

= \Nu+ md){'ip'ip) . (25) 

The fact that QGD describes confinement as well as spontaneous breaking of chiral 
symmetry is not evident from the QGD Lagrangian or the perturbative treatment 
of strong interactions described by it. It requires a non-perturbative analysis - 
lattice QGD calculations - to firmly establish the confining and symmetry breaking 
features of QGD. 

Asymptotic freedom of QGD suggests that non-perturbative effects are sup¬ 
pressed at high temperatures and hot strong-interaction matter approaches ideal 
gas behavior asymptotically. It thus is expected that non-perturbative condensates 
also disappear at high temperatures. For exact global symmetries of QGD this is 
expected to happen through a true phase transition. Basic features of this transition 
can be understood by invoking universality arguments.!!^ 

The nature of the QGD transition depends crucially on the values of the quark 
masses and the number of flavors (Nf). Fig. shows a sketch of the nature of the 
QGD transition as functions of the quark masses for a theory with two degenerate 
light (up and down) quarks with masses, mu,d = mi = mu = md^ and a heavier 
strange quark with mass, rUg, at zero baryon chemical potential. In the limit mi ^ 
oc and mg oo fermions decouple and the thermodynamics of a pure SU (3) gauge 
theory {Nf = 0) is recovered. The SU{3) gauge theory has an exact Z{3) symmetry 
and the deconflmement transition is first order. This first order region extends to 
lower quark masses and ends at a second order critical line that belongs to the 
universality class of the 3-d, Z(2) symmetric Ising modellHMl For three degenerate 
flavors Nf = 3 with small quark masses m/ = ^ 0 the chiral transition is 
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Fig. 3. A sketch of the nature of the QCD transition as functions of the two degenerate light (up 
and down) quarks with masses, = mi, and a heavier strange quark with mass, m^, at zero 

baryon chemical potential. 

known to be first order IHMl Recent lattice QCD studie^HEH with improved actions 
suggest that the extent of this first order region is quite small, i.e. limited to 
Ffii = Ffis mf^^^/270 where is the physical value of the strange quark mass. 

An additional ingredient in the discussion of the order of the transition in the 
chiral limit arises from the role of the axial anomaly. The nature of the chiral 
transition for the massless Nf = 2 theory , i.e. for mi ^ 0 and mg ^ oc, depends 
on the magnitude of the axial Ua{^) symmetry breaking. If this remains significant 
close to the transition temperature then the relevant symmetry becomes isomorphic 
to that of the 3-d 0(4) spin model and the transition is expected to be second order 
belonging to that universality class 1121££1 However, if Ua{^) symmetry breaking 
becomes negligible near the chiral transition temperature, the relevant symmetry 
becomes isomorphic to 0(2) x 0(4) and the transition be either first ordei^I^ or 
second order 1I21Z1I In the intermediate quark mass region there is no true phase 
transition, rather a crossover takes place from the hadronic to the quark-gluon 
plasma phase. 

All the first order regions are separated from the crossover region by lines of 
second order phase transitions belonging to the 3-d Z{2) universality class. The 
first order region for the Nf = 2 + 1 case, the second order Z{2) line separating 
the Nf = 2 + 1 first order and the crossover regions and the second order 0(4) line 
for the Nf = 2 case are supposed to meet at a tri-critical point characterized by a 
certain value of the strange quark mass, Although, it is well established that 

in the real world, i.e. for the physical values of the quark masses, the transition 
is a crossover pmH the location of the physical point with respect to has not 

been established and even oo cannot be ruled out. More specifically, it 
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is yet unclear whether or or < rn^p^. If 

mpys > then in the limit of m/ ^ 0 one should find a second order transition 

belonging to the 3-d 0(4) universality class, if = mp^ the tri-critical point 

is a Gaussian fixed point of the 3-dimensional model and its critical exponents 
take the mean field value^^ and if mpy^ < mp^ then in the mi ^ 0 limit one may 
cross through a second order transition belonging to the 3-d Z{2) universality class 
and then may end up in the first order transition region. 

In the following discussion we assume > mp^. 


3.1. Chiral transition in (2+1)-flavor QCD 


In the vicinity of the chiral phase transition, the free energy density may be ex¬ 
pressed as a sum of a singular and a regular part 


T 

f = -+nz 


fsing{t') P freg {T,mi,ms,il) . 


(26) 


The parameter h represents the dimensionless explicit chiral symmetry breaking 
(magnetic) field and t incorporates all the ‘thermal’ variables that do not explicitly 
break the chiral symmetry, e.g. at leading order. 


1 

to 


T - TP 

TP 




m 


and 


, 1 m; 

^ = h - 

ho rris 


(27) 


where /Xg = = jid- T° denotes the (unknown) phase transition temperature in the 

chiral limit mi ^ 0 and for fiq = 0. The scaling variables t, h are normalized by two 
unknown parameters to and ho. The other unknown parameter Kq controls how the 
chiral transition temperature T^ changes as a function of Rq. All these 4 unknown 
quantities (Tc,to, ho) are unique to QCD, similar to the low energy constants 
in the chiral Lagrangian. These parameters can be determined by analyzing the 
scaling behaviors, arising from the singular part of the free energy, of the chiral 
order parameter and susceptibilities using lattice QCD. 


3.1.1. Pseudo-critical behavior at vanishing chemical potential 

Sufficiently close to the QCD chiral transition the renormalization group invariant 
dimensionless order parameter M^, constructed out of the light quark chiral con¬ 
densate, defined in Eq. and the strange quark mass obeys the scaling relation 

Mb = = h^/^fQ^z) + regular terms , (28) 

in terms of a single scaling variable 2 : = t/h}!^^. The critical exponents [3 and 5 and 
the scaling function fop) uniquely characterize the universality class of the chiral 
phase transition. As discussed before, for continuum QCD the relevant universality 
class is expected to be the same as that of 3-d 0(4) spin models. However, the 
situation is more subtle for QCD on the lattice. For example, since for staggered 
fermions away from the continuum limit there is only one Goldstone boson in the 
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chiral limit the relevant universality class is that of the 3-d 0(2) spin model. Table[^ 
summarizes the critical exponents and other relevant quantities for the 3-d 0(2) and 
0(4) universality classes. 

Fig.j^shows results for the chiral order parameter calculated in (2H-l)-flavor 
QCpJZl] several values of the light to strange quark mass ratio mi Inis while 
keeping rUg fixed to its physical value In the continuum limit the ^ysical 

value of this ratio is = 1/27.5.1^ The left hand panel in Fig. shows 

the variation of with temperature calculated on rather coarse lattices {N^ = 4) 
for several values of mijms- In the continuum limit the smallest of these ratios 
would correspond to pion masses of about 80 MeV. Clearly, the transition sharpens 
as the light quark masses become smaller. The right hand panel of Fig. shows a 
comparison of the scaled order parameter (see Eq. 28) with the 0(2) scaling 

function. As can be seen, for small enough light quark masses, mi < mf^^^/20, the 
scaled order parameter for different values of the light quark mass collapses on a 
unique curve and compares well with the 0(2) scaling function. While these lattice 
results suggest that the chiral transition for Nf = 2 theory belongs to the 3-d 0{N) 
universality class one should bear in mind that these results are not continuum 
extrapolated and discretization effects may change the conclusion. In fact, lattice 
studies performed with the unimproved staggered action even further away from the 
continuum limit suggest that in this case the transition may even be first order.l^ 

The chiral order parameter can be used to define two susceptibilities obtained 
by either taking a derivative with respect to ‘thermal’ variable t, which give the 
mixed susceptibility Xt dMb/dt - d^fjdtdh, or with respect to the symmetry 
breaking variable, Xh ^ dMij/dh ^ d^fs/dh^. These two susceptibilities are the 
only two second derivatives of the free energy of an 0{N) symmetric theory like 
QCD that diverge at the critical poinf[^ i.e. at T = for mi = 0. Their singular 
behavior is controlled by two scaling functions, /^(^) = dfG{z)/dz, and fxi^) = 
(/g(^) ~ (^//^)/g(^))/^’ respectively. Away from the chiral limit, these two scaling 
functions have maxima at universal values of z, i.e. z = Zt and z = Zp (see Tab.j^, 


^In general, there is a third susceptibility, the specific heat, Cy, which is obtained as a second 
derivative of the free energy with respect to the ‘thermal’ variable t, Cy ~ In 0(N) 

symmetric theories, however, the relevant critical exponent a is negative. The specific heat, thus, 
does not diverge at Tc in the chiral limit. 


Table 1 . Critical expone nts a , /5, 7 , S for the 3-dimensional 
0{N) universality classesEHIl Only two of the four critical 
exponents are independent. They are related to each other 
through the scaling relations 7 = ^((5 — 1 ) and a + 2 ^ + 7 = 2 . 
The last two columns give the location of the maxi ma of the 
scaling functions /^(^) and fxi^)^ respectively EES 


N 

a 

13 

7 

(5 

Zt 

Zp 

2 

-0.017 

0.349 

1.319 

4.779 

0.46(20) 

1.56(10) 

4 

-0.213 

0.380 

1.453 

4.824 

O.TS^O) 

1.35(3) 
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Fig. 4. 0{N) scaling of the chiral order parameter Mij introduced in Eq. for (2+l)-flavor 

QCD.EHI 


which can be used to define pseudo-critical temperatures. 

For example, the scaling behavior of the chiral susceptibility for two light, de¬ 
generate flavors. 


Xm,l{T) = 2 


dmi 


(29) 


is related to that of Xh- The renormalization group invariant product of the chiral 
susceptibility Xm,i and the square of the strange quark mass is related to the scaling 
functions by. 


Vx(^) + regular terms . (30) 

At vanishing quark chemical potential the chiral susceptibility thus diverges at the 
chiral critical temperature in the chiral limit h ^ 0. The peak location of the 
chiral susceptibility will then be associated with the peak of the scaling function 
/^(z), located at z = Zp, and defines the chiral pseudo-critical temperature Tc for 
h > 0, 


1 + ^ 
zo 



+ regular terms , 


(31) 


with zq = . As illustrated in Fig. the chiral crossover temperature for 

physical QCD with rris = and j21 .h was studied in detaiP^by 

using such a scaling analysis of the chiral susceptibility, including the influence of 
regular terms. By performing scaling fits to the chiral susceptibility, extrapolating 
to the physical value of light quark mass and subsequently taking the continuum 
limit one obtains a value of 


Te = (154 ±9) MeV (32) 

for the chiral crossover temperature in (2+l)-flavor physical QCD.l^ Studies with 
other improved staggered fermion discretization schemeswhich used other 
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Fig. 5. (Left) 0(4) scaling of the chiral susceptibility and its extrapolation to the physical light 
quark mass for (2+l)-flavor QCD.^^The peak location of the chiral susceptibility defines the chiral 
crossover temperature. (Right) Continuum extrapolation of the chiral crossover temperature for 
(2+l)-flavor physical QCD, yielding Tc = (154 zb 9) MeV.^^ 

criteria, not directly related to criticality, to define a crossover temperature, and 
with chiral Domain Wall Fermion^^ also yielded compatible results for the QCD 
chiral crossover temperature. 


3.1.2. Curvature of the pseudo-eritieal line 


The scaling behavior of the mixed susceptibility, Xt^ can be utilized to determine 
the chiral phase boundary in the T-fi plane.l^ As introduced in Eq. to leading 
order the chemical potential fiq does not contribute to the explicit chiral symmetry 
breaking field h and can be treated as part of the ‘thermal’ scaling variable t. 
Thus, the derivative of the order parameter with respect to the chemical potential, 
Xm,g = is equivalent to the mixed susceptibility, Xt- The scaling 

behavior of the renormalization group invariant combination of Xni,g multiplied 
with the strange quark mass, which combines a ‘thermal’ and ‘field-like’ derivative 
of the free energy, is given by 

9k 

-(i-/3)//3(5 ^ regular terms . (33) 




2^2 


= ^h- 
to 


In the chiral limit, h ^ 0, this susceptibility diverges at for jUq = 0. For /i^ > 0, 
the chemical potential dependence of the chiral transition temperature T^{jiiq) is 
defined by the peak location, 2 : = Zt, of the scaling function and in the 

leading order in (/i^/T)^ it is given by 




1 - ^ 


^ \ 2^0 



+ regular terms . 


(34) 


As illustrated in Fig. (left), the curvature of the chiral phase transition line, 
can be determined by fitting the scaled susceptibility rUstoh^^' to the 

universal scaling function Such a scaling analysis provides a value of the 
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Fig. 6. (Left) 0{N) scaling of the susceptibility Xm,fx and determination of the curvature, Kq, of 
the chiral phase boundary in the T-p,g plane.I^ (Right) Chiral crossover temperature as function 
of baryon chemical potential, Tc(/r^), compared with the freeze-out temperature in heavy-ion 
collisions. The parametrization of the freeze-out line is taken from Ref. |86| with the freeze- 
out temperature in the limit of infinite collision energy adjusted to the crossover temperature 
at vanishing chemical potential, Tfo{\/s oo) = 154 MeV. The black band for the crossover 
line, Tc{iib)i reflects the systematic uncer tainty arising from a factor of 2 difference in current 
estimates from calculations with realHIlll and imaginar 3 EIlll fiB (see text). Data points are 
from a reweighting analysiJ^ and Taylor expansiorP^ (see text). They also are normalized to 
Tc = 154 MeV. 


curvature Kq = 0.059(6).!^ In terms of the baryon chemical potential Rb = 3/ig, 
the chiral transition temperature is then given by 


T^{rb)=T^ 


1 - 0.0066(7) 



2 

+ 0 



(35) 


This result is in very good agreement with the independent determination of the 
continuum extrapolated results of the curvature of the chiral phase boundary for 
(2+l)-flavor of physical QCDp^ using different criteria. However, analytic continu¬ 
ation from purely imaginary chemical potential yields a factor two larger values for 
the curvature EMI Fig. (right) shows the proximity of the QCD phase boundary 
to the freeze-out line in heavy-ion collisionJ^ in the T-rb plane. 


3.1.3. The critical end point 

Information on the location of a critical end point (CEP) in the QCD phase dia¬ 
gram (Fig. from lattice QCD is ambiguous. Although considerable progress in 
developing algorithms that would allow direct calculations in lattice QCD ai rb >0 


have been made recently (see Ref. 36 for a recent review) these techniques are not 
yet applicable to realistic QCD parameter values. 

The first calculations that provided hints for the existence of a critical poinfP^ 
used a reweighting technique that is applicable only on rather small lattices. It 
had been performed on coarse lattices with temporal extent = 4 using only the 
naive l-link staggered fermion action. Taste violation effects thus are large and the 
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reweighting is known to fail when used too far away from the parameters actually 
used in the calculations. It has been argued that this ‘overlap problem’ may have 
led to the spurious identification of a signal for a critical pointCalculations 
performed with an imaginary chemical potential do not find any evidence for the 
existence of a critical point .1^ However, also these calculations have been performed 
on coarse lattices and utilized the naive staggered 1-link action. 

The conclusions drawn in Ref. are based on an analysis of the /r^-dependence 
of the location of the boundary line that separates the first order transition region 
shown in Fig. at small quark masses from the crossover region. These calcula¬ 
tions, performed with staggered fermions using an imaginary chemical potential,!^ 
suggest that this boundary moves to smaller quark masses with increasing chemi¬ 
cal potential. This disfavors the existence of a critical point connected to the Z(2) 
boundary line at vanishing chemical potential. A similar analysis, performed with 
clover-improved Wilson fermions, comes to the opposite conclusioiP^ and thus does 
favor the existence of a critical point. The apparent differences between calcula¬ 
tions performed with staggered fermions and Wilson fermions on coarse lattices 
underscores that better control over systematic errors is needed before a definite 
conclusion on the existence of a critical point can be drawn. 

Analyzing the convergence properties of the Taylor series for the pressure 
at non-zero chemical potentiaEnEH (Eq. 18) does, in principle, allow to relate 
the location of the CEP to the radius of convergence of this series. Although 
this approach has been used closer to the continuum limit {N^ = than 

the studies discussed above, the current estimates still are based on calcula¬ 
tions with the naive staggered l-link action. They still suffer from large taste 
symmetry violations. These calculations yield for the location of the critical 
poinfP^ (T^/Tc,yuf/T^) = (0.94(1), 1.68(5)which is significantly lower than 
the reweighting result {T^/Tc^ = (0.99(1), 2.2(2)).!^ We show these two 

estimates for the location of the CEP in Eig. (right). The systematic errors of 
these estimates are at present difficult to estimate. 

Eurther information on the location of a critical point comes from the calcu¬ 
lations of (i) the equation of state at /i^ > 0 performed with improved staggered 
fermions, which we will discuss in Section]^ as well as (ii) the analysis of freeze-out 
parameters based on Taylor expansions of cumulants of charge fluctuations, which 
are also performed with improved staggered actions and which we will discuss in 
Section]^ both suggest that a CEP located at /r^/T < 2 is unlikely. 


3.2. Deconfining aspect of QCD 

Pure SU{Nc) Yang-Mills theory without quarks possess an exact global Z{Nc) cen¬ 
ter symmetry, which gets spontaneously broken in the high temperature deconfining 
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phaseThe thermal expectation value of the renormalized Polyakov loop 


iren(T) - 6 ( '^ 11 U(xo,X),6 ) ’ 

^ ^ \ xo = l / 


N, 


(36) 


is not invariant under the Z{Nc) center symmetry. The Polyakov loop can be in¬ 
terpreted as the free energy difference, Foo(T), of a thermal system with and with¬ 
out an infinitely heavy static quark anti-quark pair separated by infinite distance, 
^ren(T) = exp(—Foo(T)/2T). The renormalization constant c{g‘^) can be fixecP^by 
demanding that in the short distance limit the heavy quark free energy coincides, 
up to a trivial additive constant, with the Coulombic short distance behavior of the 
zero temperature heavy quark potential defined in Eq. 23 In the confined phase 
Foo(T) = oo, as a static quark and anti-quark pair cannot be separated by infinite 
distance in this phase. Thus in the Z{Nc) symmetric confined phase Lj.q^{T) = 0. 
On the other hand, in the spontaneously Z{Nc) broken deconfined phase a static 
quark anti-quark pair can be separated by infinite distance due to the presence of 
color screening and Lj.q^{T) ^ 0. Thus, for a pure SU{Nc) gauge theory, i.e. in the 
limit m/ ^ oc, the Polyakov loop serves as the order parameter for the deconfine¬ 
ment transition.!^ However, since the mere presence of quarks explicitly breaks the 
Z{Nc) center symmetry, the Polyakov loop does not serve as an order parameter 
for QCD with realistic light quarks. Furthermore, since the Polyakov loop is not 
related to a derivative of the QCD partition function with respect to the thermal 
or the symmetry breaking field, its change or fluctuation across the QCD transition 
may not capture the true singularities of the QCD partition function in any limit. 
Thus, a deconfining temperature defined from the change of the Polyakov loop may 
not reflect the pseudo-critical properties of QCD with realistic light quark masses. 
As shown in Fig. 0(1 eft), the change of the Polyakov loop within the chiral crossover 
region, Tc = 154(9) MeV, is rather gradual and smooth. 
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Fig. 7. (Left) The renormalized Polyakov loop in (2+l)-flavor QCD.I^ (Right) Appearance of the 
fractionally charged degrees of freedom in the chiral crossover region Tc = 154(9) MeV (shaded 
region) for the light as well as the strange quark. The black points show results obtained using 
the stout actiorff^ and the other points have been obtained using the HISQ action.!^ 
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Deconfinement is associated with the liberation of degrees of freedom, mani¬ 
fested by a rapid rise in the bulk thermodynamic observables such as the pres¬ 
sure, energy density etc. Among such bulk thermodynamic observables fluctuations 
and correlations of different conserved quantum numbers, e.g. the baryon number 
(5), electric charge (Q) and strangeness (S'), directly probe the liberation of the 
quark degrees of freedom and the consequent appearance of fractionally charged 
quantum number carriers. For example, both for an uncorrelated hadron gas and 
in a free quark gas the difference of second to fourth order cumulants of baryon 
number fluctuations probes the baryon number of degrees of freedom, specifically 
X 2 — Xa ^ For hadronic degrees of freedom carrying baryon number 

B = ±.1 one thus has X 2 ~ Xa — O 5 t)ut for quark-like degrees of freedom with 
B = ±1/3 one finds X 2 ~Xa 7^ 0- Similarly, for the strange quark sector higher or¬ 
der baryon-strangeness correlation such as X 31 ~X 11 ^ {B^ — B)S is non-vanishing 
if the strangeness carrying degrees of freedom are associated with quark-like baryon 
number B = ±1/3 and vanishes when the strangeness is carried by hadronic degrees 
of freedom with 5 = ±1. Thus, such combinations of higher order fluctuations and 
correlations of quantum numbers are sensitive probes of appearance of fraction¬ 
ally charged degrees of freedom in a deconfined medium EUUMI xhe right hand 
panel of Fig. shows lattice QCD results for X 2 ~ Xa Xsf ~ Xif from the 
BNL-Bielefeld collaboratiorP^ as well as the Budapest-Wuppertal collaboration.^^ 
Clearly, both quantities start deviating strongly from zero in the chiral crossover 
region, Tc = 154(9) MeV, indicating that fractionally charged degrees of freedom 
start to appear at these temperatures and onset of deconfinement takes place for 
the light as well as the strange quarks. Note, however, that also these fluctuation 
observables are not order parameters in the strict sense. 

3.3. Axial symmetry of QCD at high temperature 

Although, the axial Ua{^) symmetry is not an exact symmetry of QCD, as men¬ 
tioned before, the magnitude of its breaking is expected to influence the order of 
the chiral phase transition. Hence, knowledge of the temperature dependence of the 
Ua(X) breaking is essential for a comprehensive understanding of the QCD chiral 
transition. The axial Ua{^) symmetry of the massless QCD Lagrangian is broken 
due to quantum fluctuations, resulting in non-conservation of axial currenfE P^F^^ I 
as well as explicit breaking of the global Ua{^) symmetry induced by topologi¬ 
cally non-trivial gauge field configurations, such as the instantons.^^^ Due to color 
screening the instanton density gets suppressed as the temperature increase^^^ and 
in the T ^ 00 limit the Ua{^) symmetry becomes exact. 

Since Ua(X) exact symmetry of QCD one cannot define an order pa¬ 

rameter associated with its breaking. However, for two massless flavors, the pion 
(tt) and the iso-vector scalar S (ai) mesons transform into each other under a Ba(1) 
rotation. Since the presence of an exact Ba(1) will render these meson states degen¬ 
erate, the difference of the integrated two-point correlation functions of pion and S 
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meson, 

X 7 r-Xs = J dx (7r+(a;)7r“(0)) - ((5+(a;)(5"(0)> ^ , (37) 

will also vanish in the limit of exact Ua{^) symmetry. Corrections due to small non¬ 
vanishing light quark masses will only contribute at Thus, this quantity 

can be taken as a measure of the Ua{'^) breaking .CSH 



Fig. 8 . The Ua{^) breaking measure Xvr — XS renormalized in the MS scheme as function of 
temperature. Results are from a calculation with Domain Wall fermions.I^ The difference is non¬ 
vanishing for T > Tc = 154(9) MeV. It becomes independent of quark mass for T > 168 MeV and 
rises with decreasing quark mass at lower temperature. In fact, it will diverge in the chiral limit 
for T < Tc. 


This measure of Ua{^) breaking was studied in detail using improved p4 stag¬ 
gered fermions ,11^ and it was found that Xn — Xs remains non-vanishing for 
T < 1.2Tc. However, for staggered fermions the issue of axial anomaly is quite 
subtle and the correct anomaly may only emerge in the continuum limit E^ZE^ll 
On the other hand, emergence of the axial anomaly is more straightforward for 
the chiral DWF formulation.l^^ For DWF, axial symmetry is broken by the same 
topologically non-trivial configurations as in the continuum. Lattice artifacts ap¬ 
pear only at due to the explicit chiral symmetry breaking residual rnass,!^ 

rures, arising from finiteness of the fifth dimension {Lg < oo). Thus, the DWF 
action turns out to be a natural candidate for investigation of the temperature de¬ 
pendence of Ua{^) breaking in QCD. The temperature dependence of the Ua{^) 
breaking measure, X'k —Xs^ has been extensively studied using the DWF formalism 
for several volumes as well as quark masses EHMEMl As depicted in Fig. cal¬ 
culations with DWF clearly show that X'k — X5 does not vanish around the chiral 
crossover temperature and remains non-vanishing for 165 MeV ^ T < 195 MeV, 
independent of the light quark masses. These results indicate that Ua{X) symmetry 
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may remain broken at these temperatures even in the chiral limit. 



A/T 


\|/+V|/ V|/^Y5V 



Fig. 9. (Top) Accumulation of near-zero eigenmodes of the Overlap-Dirac fermion matrix at 
T = 1.5Tc.^^ (Bottom) Spatial profile of a typical near-zero mode (left) and the spatial profile of 
the chirality of the same near-zero mode (right) at T = 1.5Tc.^^ 


Since Ua{^) breaking arises due to topology of the gauge fields, it is intimately 
related to the infrared modes of the Dirac fermions. In the limit of infinite volume, 
both the chiral order parameter, ('^'0)^ and the Ua{^) breaking measure, X'k — X5^ 
can be represented in terms of the eigenvalue density, p(A), of the Dirac fermions 



2mi p{X) 
+ mf 


and 


Xtt - X(5 = 



p(A) 

(A^ + mf)^ 


(38) 


In the chiral symmetric phase and in the chiral limit must vanish, but Xtt ~X(5 

can remain non-zero. Identification of the infrared fermionic modes, i.e. the form 
of p(A), that give rise to such a phenomenon naturally leads to the underlying non- 
perturbative mechanism of axial symmetry breaking. Studies with DWFl^Hil^ as 
well as with overlap fermions ,1^ possessing even better chiral properties and an exact 
index theorem, suggest that an accumulation of near-zero eigenmodes of the form 
p(A) ^ m‘f5{X) may largely account for the observed Ua{^) breaking in Xtt — X6 at 
high temperature. Such accumulation of the near-zero modes are depicted in Fig. 1^ 
(top). More detailed studied of the space-time profiles, as illustrated in Fig. 
(bottom), localization properties and distributions of these near-zero modes indicate 
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that their behavior is consistent with underlying presence of a dilute instanton gas— 
a gas of widely separated, weakly interacting, small instantons and anti-instantons. 
This suggests that even at temperatures T ^ 1.5Tc weakly interacting instanton 
anti-instanton pairs are largely responsible for the Ua{^) breaking. 

The analysis of symmetry breaking close to the chiral limit of 2 or (2+1)- 

flavor QCD, however, is far from being settled. Detailed systematic studies of 
the quark mass and volume dependence as well as the cut-off dependence are still 
missing. Moreover, results obtained with chiral fermions are still controversial. 
Calculations performed with dynamical overlap fermionand so-called optimal 
domain wall fermionJ^^ at present suggest that Ua{^) does get restored already 
at Tc, which would still allow for the occurrence of a first order transition when 
approaching the chiral limit. 

4. Bulk thermodynamics 

4.1. The QCD equation of state at vanishing chemical potential 

Lattice QCD calculations of the equation of state, or more general, of basic bulk 
thermodynamic observables like the pressure (P), energy density (e) or entropy 
density ( 5 ), have reached now a precision where continuum extrapolated results 
for physical light and strange quark masses can be obtained. This also allows 
detailed comparisons with perturbative calculations at high temperature and model 
calculations at low temperature. 

The procedure to calculate bulk thermodynamic observables in lattice QCD 
is well established. The starting point for these calculations is the evaluation of 
the trace anomaly, 0^^(T) = e — 3P. The trace anomaly is directly obtained as a 
temperature derivative of P/T^ which also is related to the QCD partition function. 


Z{T,V) 


e-3P _^dPlT^ 
T4 “ dT 


(39) 


with P/T^ = \imv^oo{VT^)~^ lnZ(T, V). As the trace anomaly arises as a deriva¬ 
tive of the logarithm of the QCD partition function with respect to temperature, 
it can be expressed in terms of expectation values of rather simple observables, e.g. 
a combination of the gauge action and light and strange quark chiral condensates. 
Nonetheless, the calculations become rather time consuming because a subtrac¬ 
tion of the zero temperature contributions is needed to eliminate divergent vacuum 
contribution^ 

As can be seen from Eq. up to an integration constant at a temperature Tq 
the pressure can be obtained from the trace anomaly. 



( 40 ) 


‘^For more details we refer, for instance, to Ref. 112 
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Fig. 10. (Left) Comparison of th e tra ce anomaly (e — 3P )/T^, pressure and entropy density 
calculated with the HISQ (coloredand stout (greyj^^ discretization schemes for staggered 
fermions. (Right) Continuum extra polat ed results for pressure, energy density and entropy den¬ 
sity obtained with the HISQ action.^^ Solid lines on the low temperature side correspond to 
results obtained from hadron resonance gas (HRG) model calculations. The dashed line at high 
temperatures shows the result for a non-interacting quark-gluon gas. 


This allows to reconstruct the energy density as well as the entropy density s/T^ = 
(e + P)/Tt 

The determination of thermodynamic quantities in QCD is a parameter free 
calculation. All input parameters needed in the calculation, e.g. the quark masses 
{rriu = rrid^ rris) and the relation between the lattice cut-off, a, and the bare gauge 
coupling, /3 = 6/^^, are determined through calculations at zero temperature. Like¬ 
wise, there is only a single independent thermodynamic observable that is calculated 
in a lattice QCD calculation, for instance the trace anomaly, 0^^(T). All other bulk 
thermodynamic observables are obtained from through standard thermo¬ 

dynamic relations. In Fig. (left) we show recent results for the trace anomaly 
of (2+l)-fiavor QCLEMHI obtained with two different discretization schemes by 
two different groups. The results are extrapolated to the continuum limit and are 
obtained with a strange quark mass tuned to its physical value and light quark 
masses that differ slightly {rus/mi = 27l^^and 2(jni)^ 

The right hand panel in this 
figure shows results for the pressure, energy density and entropy density obtained 
from the trace anomaly by using Eqs. ^ and[^ 

Also shown in Fig. 10 are results obtained from a hadron resonance gas (HRG) 
model calculation of bulk thermodynamics. As can be seen this describes the QCD 
equation of state quite well also in the transition region, although it may be noted 
that the HRG calculations yield results for all observables that are at the lower error 
band of the current QCD results. It has been speculated that this may indicate 
contributions from additional, experimentally not yet observed resonances which 
could contribute to the thermodynamics.^!^ Indeed evidence for the contribution 
of a large number of strange baryons has recently been found in lattice QCD calcu¬ 
lations of conserved charge fluctuation^!!^ (see also the discussion in Section and 

0 - 
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Fig. 11. The critical energy density Cc in (2+l)-flavor QCD. The band gives the continuum 
extrapolated result for the energy density taken from Ref. EH] The HRG curve is based on all 
resonance with mass less than 2.5 GeV listed by the Particle Data Group.^^ 


In Fig. [H we show the energy density in the low temperature region. The 
box highlights the transition region characterized by the crossover temperature, 
Tc = (154 ±9) MeV. From this we deduce the energy density in the crossover region, 
ec = (0.34 ± 0.16) GeV/fm^. This is a rather small value for the energy density 
needed to convert ordinary hadronic matter into a medium made up from quarks 
and gluons. It may be compared to the energy density of ordinary nuclear matter, 
^nuclear matter ^ q GcV/fm^ or the energy density inside a nucleon, 

0.45 GeV/fm^, assuming the radius of nucleon Rn — 0.8 fm. In fact, ec is close to 
the energy density reached in the dense packing limit of nucleons with radius 

The simple bulk thermodynamic observables like pressure and energy density 
give the impression that an HRG model also provides a good description of the 
equation of state for temperatures above the crossover region, i.e. for T > 165 MeV. 
However, as discussed in the previous section the analysis of conserved charge fluc¬ 
tuations clearly shows that at temperatures T > 160 MeV thermodynamics can 
no longer be described in terms of hadronic degrees of freedom (see for instance 
Fig. 0 (right) and the discussion in Section]^. This also becomes evident in sec¬ 
ond order derivatives of the QCD partition function with respect to temperature. 
The speed of sound, = dp/de, is related to the inverse of the specific heat, 
Cv = de/dT, 


Cy 

T3 


dp dp/dT s 
d^ ^ de/dT ^ ^ ’ 


T _ G e d{e/y) 

dT y V T4 dT 



(41) 

(42) 


Both quantities are shown in Fig. The specific heat does not develop a pro¬ 
nounced peak in the transition region as one could have expected from pseudo- 
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Fig. 12. The velocity of sound in (2+l)-flavor QCD (left) and the specific heat Cy jT^ together 
with the two components (see Eg. |42| ) contributing to it (right). Solid black lines in the low and 
high temperature regions show the corresponding hadron resonance gas (HRG) and non-interacting 
quark-gluon gas results, respectively. 

critical behavior of energy density fluctuations close to a critical point. This may 
be understoocjnil from the temperature dependence of the two terms contribut¬ 
ing to Cy/T^- The dominant singular contribution arises from the temperature 
derivative of e/T^, which has a peak. This, however, is overwhelmed by the large 
energy density contribution at high temperature which reflects the liberation of 
many partonic degrees of freedom. Furthermore, even in the chiral limit, where 
QCD is expected to have a second order phase transition belonging to the uni¬ 
versality class of 3-d, 0(4) symmetric spin models, the specific heat will not di¬ 
verge as the relevant critical exponent a —0.2 that controls its singular behavior, 
CvjT^ ^ (|T — Tc\/Tc)~^ + const., is negative for this universality class (see Ta¬ 
ble [^. The speed of sound will therefore stay non-zero at Tc also in the chiral limit. 


4.2. The QCD equation of state at non-vanishing chemical potential 

In Section we have discussed lattice QCD results on the dependence of the QCD 
crossover temperature on the baryon chemical potential and its relation to the 
freeze-out temperatures determined in heavy ion experiments. These experiments, 
in particular the beam energy scan program performed at RHIC, will probe prop¬ 
erties of strong-interaction matter at non-vanishing baryon chemical potential in 
the temperature range 0.9 < T/Tc ^ 2 and 0 < ^b/T < 3, with Tc denoting 
the crossover temperature at = 0. For the hydrodynamic modeling of matter 
in this (T^iib) regime it thus is of importance to also know the equation of state 
at non-vanishing iib/T. As direct numerical calculations at non-zero fiB are not 
yet possible, a viable approach is to analyze the equation of state using a Taylor 
expansion in terms of chemical potential^HESl as given in Eq. 


18 


In this way 

some results for the EoS at non-zero baryon chemical potential have already been 
obtained on coarse latticesE31IEIlIIl2l Continuum extrapolated results for Taylor 
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expansion coefficients of the pressure at have been obtained EMUIl These 

have been used to construct the EoS at which implements the strangeness 

neutrality constraints^ {{^s) = 0) and the electric charge to baryon number rela¬ 
tion {{uq) = 0 . 4 ( 77 , 5 ) ) suitable for conditions met in heavy ion collisions. 




Fig. 13. Expansion coefficients of the pressur e at non-zero baryon chemical potential. The left 
hand figure shows the leading order correctiorS^ and the right hand figure shows the relative 
contribution of the next to leading order correction. The continuum extrapolated result obtained 
with the stout action is taken from Ref. CMl 


We will discuss here only the case rq = rs = 0 in some detail in order to 
illustrate the relative importance of higher order corrections in different temperature 
and rb regions. The Taylor series for the pressure is given by, 


P{T,fiB)-P{T,0) 

2^4 




1 X?{T) 

12 xiiT) 



+ . 


(43) 


The leading order correction to the pressure thus is proportional to the quadratic 
fluctuations of net baryon number. The next to leading order corrections are pro¬ 
portional to the quartic fluctuations. When written in the form given in Eq. 
it becomes easy to identify the relative importance of leading and next to leading 
order corrections. In Eig. 13 we show xf (T) (left) and xf (^)/xf (^) (right). With 


increasing temperature the 0{rq) correction rapidly looses importance relative to 
the leading 0{r\) term. 

In Eig. we show preliminary results for the /i 5 -dependent contribution to 
the total pressure evaluated for different values of rb/T and taking into account 
corrections up to 0{{rb/TY)^^ These results suggests that an 0 {{rb /TY) Taylor 
expansion of the pressure (and energy density) is well controlled for all values of the 
chemical potential below rb/T = 2. This covers a wide range of the QCD phase 
diagram accessible in the beam energy scan (BES) at RHIC, i.e. the region of beam 
energies x/s > 20 GeV. 
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Fig. 14. (Le ft) The /x^-dependent part of the pressure at 0{{faB/T)^) (black) and 0{{iib 
(colored) The latter is shown only in the temperature regime where neglected corrections at 
/T)^) are estimated to contribute less than 10%. (Right) When combined with the jiB = 0 
contribution for t he p ressure shown in Fig. |10| these neglected terms contribute less than 3% to 
the total pressureThe grey band gives the uncertainty on_P(T, 0)/T^ and the central line in 
the band is the parametrization of P(T, 0)/r‘^ given in Ref. |ll4] 


4 . 3 . Perturbation theory, hadron resonance gas and the strongly 
interacting liquid 


In our discussion of bulk thermodynamic observables we have compared with HRG 
model calculations. We have seen that the HRG does provide a rather good de¬ 
scription of bulk thermodynamics below and even in the crossover region. However, 
the HRG does depend on the resonance spectrum used in the calculation. This is of 
some significance when one looks at more selective observables like, e.g., fluctuations 
and correlations of strange baryons (see Sections and [^. Spectrum independent 
HRG results are, for instance, related to fluctuations of net baryon number. The ra¬ 
tio of fourth to second order net-baryon number cumulants shown in Fig, (right) 
is unity in any HRG model irrespective of details of the spectrum. This figure thus 
shows that an HRG model calculation can be expected to provide a good descrip¬ 
tion of QGD thermodynamics on a 10% level only at temperatures less than about 
T - (140 - 145) MeV. 

At asymptotically high temperatures hard thermal loop, perturbative calcula¬ 
tion^^ or dimensionally reduced QGD (EQGD)P^ provide a well-established frame¬ 
work for the analysis of bulk thermodynamics and fluctuations of conserved charges. 
We show in Fig. 15 a comparison of the trace anomaly (left) and the pressure 
(right) calculated in lattice QGD with HTL and EQGD calculations. These cal¬ 
culations seem to give a good description of bulk thermodynamics at tempera¬ 
tures T > 400 MeV. However, even in this regime there remain differences between 


HTL and EQGD calculations that are of the order of 10%. Moreover, within the 
HTL calculations there seems to be no unique choice for the renormalization scale 
/i, that would allow to match all observables simultaneously. Similar conclusions 
can be drawn from the HTL analysis of net-baryon number fluctuations shown in 
Eig. (right). Although observables that are dominated at high temperature by 
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Fig. 15. Comparison of the (2+l)-flavor calculatiorP^^ of the trace anomaly (left) and pressure 
(right) with HTL and EQCD (dashed line) calculations. The black line corresponds to the HTL 
calculatiorPS with renormalization scale fi = 2 ttT. Note that this solid line would move up for the 
trace anomaly and move down for the pressure if the scale fi in HTL is reduced. 




Fig. 16. The ratio of quartic and quadratic net-baryon number fluctuations versus temperature. 
The left hand panel shows temperature ranges in which HRG and resummed perturbative calcu¬ 
lations, respectively, provide good approximations to lattice QCD results. The right hand panel 
shows the result from a HTL-resummed calculations.!^ 

quark rather than gluon contributions seem to approach perturbative behavior ear¬ 
lier, it still is evident that agreement with lattice QCD results on the 10% level only 
is possible for T > (250 — 300) MeV. In general the temperature range Tc <T <2Tc 
is highly non-perturbative and obviously not accessible to hadronic model calcula¬ 
tions. This is highlighted in the left hand panel of Fig. We will discuss in the 
following sections properties of strong-interaction matter in this temperature range. 


5. Fluctuations of conserved charges 

Proximity of a second order criticality, such as the 0(4) chiral phase transition 
or the QCD critical point, is universally manifested through long-range correla¬ 
tions at all length scales, resulting in increased fluctuations of the order parameter. 
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These fluctuations can be quantified through the Gaussian (variance) as well as 
non-Gaussian (skewness, kurtosis etc.) cumulants of the distribution of the order 
parameter. Higher order non-Gaussian cumulants become increasingly sensitive to 
proximity of a critical poinfP^ which is reflected in their grows with higher powers 
of the correlation length.Moreover, even qualitative features, such as the sign 
change and the associated non-monotonicity, of these non-Gaussian cumulants can 
encode the presence of a nearby critical region.Non-Gaussian cumulants can 
be accessed in heavy-ion experiments via the event-by-event fluctuations of vari¬ 
ous conserved charges and particle multiplicities.fl^^®^ In this vein, a major focus 
of the Beam Energy Scan (BES) progranP^ at the RHIC is measurements of the 
event-by-event fluctuations of particle multiplicities and conserved charges 

Among several conserved charges, the net electric charge is of special interest. 
Gumulants of net electric charge fluctuations can be measured in experiment^MEHl 
and are calculable in lattice QCDEMHII Gumulants of net electric charge fluctua¬ 
tions and correlations with other conserved charges are well defined in lattice QGD 
at vanishing chemical potentials where standard lattice QGD techniques can be used 
to compute them. Eurthermore, Taylor expansions in powers of around /iB = 0, 
can be employed to obtain generalized susceptibilities for >0, 

OO 1 , 

= . (44) 

/c=0 


These susceptibilities directly relate to cumulants of the net electric charge fluctu¬ 
ations. 


X? 

X? {T,^b) 
X2 {T,^ib) 
X? {T,^ib) 


1 

vF 

1 

vF 

1 

vF 

1 

vF 


(Nq) , 

{mQf), 
{mQf), 
{mQ)^) - 




(45) 


where Nq is the net (positive minus negative) charge and SNq = Nq — (Nq). 

On the other hand, through the measurements of the event-by-event distribu¬ 
tions of the net electric charge, heavy-ion experiments provide various cumulants, 
mean (Mg), variance (erg), skewness (S'g), and kurtosis (/^g), of the electric charge 
fluctuations for given beam energy 


Mq (^/i) = {Nq) , 


Sq{V^) 



= (rngf) 


kq (Vs) 



(46) 


Thus, the charge susceptibilities obtained from lattice QGD calculations and the 
cumulants measured in the heavy-ion experiments are directly related to each other 
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through the appropriate volume-independent ratioJ^^ 

Mq {^/s) ^ Xi {T, Mb) ^ 
X?(T,/xb) 

(aA) (aA) ^ X3 {T, Mb) ^ jjQ 
Mq{V~s) ~ Xi{T,I^b)~ 


(47a) 

(47b) 


As in many other cases, one still has to worry about directly confronting thermody¬ 
namic QCD calculations performed in equilibrium and in the thermodynamic limit 
to the rapidly expanding dense matter created in heavy ion collisions in a small 
volume. In the analysis of higher order cumulants additional complications may 
arise from more complicated efficiency corrections and the limited experimental 
acceptance windows IIMIMSJ 

In heavy-ion collision experiments the only tunable parameter is the beam en¬ 
ergy, x/s. However, to gain access to the information regarding the QCD phase 
diagram this only tunable parameter needs to be related to the thermodynamic 
variables, e.g. temperature and baryon chemical potential. Traditionally, this 
{T^rb) mapping has been done by relying on the statistical hadronization 
model based analysis.l^ Recent advances in heavy-ion experiments as well as in 
lattice QCD calculations have placed us in a unique situation where, for the first 
time, this mapping can now be obtained through direct comparisons between the 
experimental results and rigorous (lattice) QCD calculations. Recently, it has been 
showrUHlli^ that by directly comparing lattice QCD calculations for (Eq. 47b) 
and (Efi- 47a) with their corresponding cumulant ratios measured in heavy-ion 
experiments it is possible to extract the thermal parameters, namely the freeze-out 
temperature, and the freeze-out baryon chemical potential The feasibility 
of such a procedure has been demonstrated in Refs. 123|l47|l'48 Fig. [Tf31 

illus¬ 
trates a recent example of such a comparison and subsequent determination of the 
freeze-out parameters. 

As can be seen from Fig. due to large errors on the experimental results 
for (5'Qcrg)/Mg, at present, only an upper limit on the freeze-out temperature 
can be determined using the method described above. Thus, a complementary 
procedure for the determination of relying on a separate observable that can 
be extracted both from heavy-ion experiments and lattice QCD calculations, is 
certainly welcome. Recently, such a complementary procedure for the determination 
of T4 has been proposed in Ref. 116 This procedure takes advantage of the fact 


that the initially colliding nuclei in heavy-ion collisions are free of net strangeness. 
Thus, the conservation of strangeness under strong interaction ensures that the 
QGP medium created during the collisions of these heavy-ions is also strangeness 
neutral. 

By Taylor expanding the net strangeness density, {ns) (rb^Rs)^ in Rb and rs 
and subsequently imposing the strangeness neutrality condition, (ns) {rb^Rs) = 0, 
for a homogeneous thermal medium the strangeness chemical potential, rs, can be 
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Fig. 17. (Left) A comparison between the lattice QCD result^^for and the STAR datcE^ 
for {Sq (7^)/Mq at = 62.4 GeV. The overlap of the experimental results with the lattice 
QCD calculations provides an upper bound on the freeze-out temperature < 155 MeV. (Right) 
Lattice QCD resultfl^for R ^2 ^ function oi /a b compared with the STAR datcff^for MqjaQ 

in the temperature range = 150(5) MeV. The overlap regions of the experimentally measured 
results with the lattice QCD calculations provide estimates for the freeze-out chemical potential 
for a given ^/s. The arrows indicate the corresponding values of ob tained from the statistical 
hadronization model fits to experimentally measured hadron yields!^ 



obtained 


Rs 

Rb 


si(T) + S3(T)(^)' + 0 



(48) 


The coefficients 5 i, 53 , etc. consist of various generalized baryon, charge and 
strangeness susceptibilities defined at vanishing chemical potentials and can be cal¬ 
culated through standard lattice QCD computations at zero chemical potentials.!!^ 
Fig. [^(left) shows the leading order contribution to Rs/rb: be. si{T). A compar¬ 
ison of the lattice result with the predictions from the hadron resonance gas model 
reveal that the inclusion of only experimentally observed hadrons, as listed by the 
Particle Data Groupfails to reproduce the lattice results around the crossover 
region. Note that, while Rs/rb is unique in QCD, for a hadron gas it depends on 
the relative abundances of the open strange baryons and mesons. For fixed T and 
/i^, a strangeness neutral hadron gas having a larger relative abundance of strange 
baryons over open strange mesons naturally leads to a larger value of rs- Aston¬ 
ishingly, the inclusion of additional, unobserved strange hadrons predicted within 
the quark modell^^llHI provides a much better agreement with lattice results, hint¬ 
ing that these additional hadrons become thermodynamically relevant close to the 
crossover temperature.!!!^ As will be discussed in Sectionj^, other lattice thermody¬ 
namics studies also indicate that additional, unobserved charm hadrons do become 
thermodynamically relevant close to the QCD crossover.!!^ 

On the other hand, the experimentally measured yields of the strangeness S 
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Fig. 18. (Left) Lattice QCD result J^^for \is jT b at the leading order, i.e. si(T') (see Eq. |48| . The 
dotted line (PDG-HRG) shows the results of a hadron resonance gas model containing only hadrons 
listed by the Particle Data Group.^^ The solid line (QM-HRG) depicts the result for a had ron gas 
when additional, experimentally yet unobserved, quark model predicted strange hadroni^2lISD 
are included. The shaded region indicates the chiral crossover region Tc = 154(9) MeV. (Right) 
A comparison between experimentally extrac ted v alues of /i^/T^) (filled symbols) and 

lattice QCD results for fis/T b (shaded bands).^^^ The lattice QCD results are shown for iibIT = 
yJ^jT^. The temperature range where lattice QCD results match with provide the values 

of T-^, i.e. T^ = 155(5) MeV and 145(2) MeV for = 39 GeV and 17.3 GeV, respectively. 


anti-baryon to baryon ratios, at the freeze-out are determined by the thermal 
freeze-out parameters (T-^,/i;^, 


RH{\fs) = exp 


Tf 



(49) 


By fitting the experimentally measured values of Rs and Rq, corresponding 
to |5'| = 1,2 and 3, the values of Rs/i^b ‘observed’ in a heavy-ion 

experiment at a given can easily be extracted. Matching these experimentally 
extracted values of r^/i^b with the lattice QCD results for rs/rb as a function of 
temperature, one can determine the freeze-out temperature T-^. Fig. 18 (right) il¬ 


lustrates this procedure. Once again, the inclusion of additional unobserved strange 
hadrons in the hadron resonance gas model (QM-HRG) leads to very similar val¬ 
ues of the freeze-out temperatures as obtained using the lattice data. However, 
including only the hadrons listed by the Particle Data Groupl^ (PDG-HRG) yields 
freeze-out temperatures that are 5-8 MeV larger. 


6. Transport properties 

As we have seen in previous sections the analysis of bulk thermodynamics and charge 
fluctuations provides plenty of evidence that thermodynamics of strong-interaction 
matter above the crossover transition temperature, Tc, and up to temperatures 
of about (1.5 — 2)Tc is highly non-perturbative. In this temperature range the 
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velocity of sound is significantly smaller than that in a non-interacting quark gluon 
gas suggesting also significant modifications of its transport properties. One of 
the striking results of heavy-ion experiments at RHIC and LHC is that strong- 
interaction matter in this temperature range shows features of a strongly coupled 
medium with a small shear viscosity to entropy density ratio and small diffusion 
coefficients. This led to the notion of strong-interaction matter close to the crossover 
transition being a ‘nearly perfect liquid’. 

Transport properties of strong-interaction matter can be explored in lattice QCD 
calculations through the analysis of current-current correlation functions, e.g. for 
some current J (r, x ), 

/ pOO J 

d3xe*PV^“(0,0)J^^(r,x)) = y^ (50) 

where the integration kernel is 


K{uj,r) 


cosh(u;(r — 1/2T)) 
sinh(cc;/2T) 


(51) 


and (a, P are appropriately chosen (scalar, vector, tensor) indices specifying the 
current. 

In Euclidean space-time these correlation functions have a spectral representa¬ 
tion. The low frequency and low momentum structure of the spectral functions 
Pj^(co’,p) then gives access to transport properties of strong-interaction matter 
while at higher frequencies the spectral functions provide information on in-medium 
properties of bound states. We will discuss the latter in the Section Extracting 
the spectral function from a Euclidean correlation function is, however, difficult. 
This is due to the fact that Euclidean correlation functions are calculated in finite 
temperature lattice QCD only at a few discrete Euclidean time separations, e.g. 
0 < T < Nt It requires the application of inversion methods like the Maximum 
Entropy Method (MEM]P^ or a modeling of the structure of the spectral function 
at low frequencies. 

According to the Kubo relations, which express the response of a medium to 
small perturbations, transport coefficients can be determined from spectral func¬ 
tions calculated in equilibrium thermodynamics. They are given by the slope of the 
corresponding spectral functions at vanishing energy. Eor instance, using in Eq. 
the electromagnetic current J^(t,x) = Jy(t,x) = ' 0 (t,x) 7 ^' 0 (t,x) for a fermionic 
charge carrier (quarks) with charge e (e.g. Qd = —e/3 for a down quark) one can 
extract the electrical conductivity a which quantifies the response of the QGP to 
small perturbations induced by an electromagnetic field. 


^ lim lim V 

6 p^O 


UJ 


(52) 


"^Note that Eq. |5l| is symmetric around r = 1/2T and as a consequence only half the number 
of data points of the Euclidean temporal correlation function, i.e. Arr/2, provide independent 
information to extract the spectral function. 
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Here Cem is the sum of the square of the elementary charges of the quark flavor 
/, Cem = Q^f‘ The electrical conductivity and the vector spectral function are 
also related to the emission rate of soft thermal photons, 


, dR^ 
hm uj 
uj^o d'^p 


= ^ (^{T) T aem 


(53) 


and the contribution of quark-antiquark annihilation to the thermal dilepton rate, 


dVk 


bai 


1 


/9y(u;,p,T), 


(54) 


dcjd^p 547r^ — 1) 

where aem is the electromagnetic fine structure constant. 

Using the vector current Jy(r, x) for heavy quarks one obtains from Eq. 50 the 
corresponding spectral function of heavy quarkonium in the vector channel, py, 
which allows to determine the heavy quark diffusion coefficient 

D=S-lim lim V ( 55 ) 

6^2 <-^^0 p^O ^^ ' • 


1=1 


UJ 


Here X 2 heavy quark number susceptibility which is defined through the 

zeroth component of the temporal meson correlator in the vector channel. E.g., for 
charm quarks this is just the net charm number susceptibility introduced in Eq. 
Similarly, in the light quark sector, the electric charge diffusion coefficient Dq is 
defined as the ratio of a to electric charge susceptibility X 2 defined in Eq. 

Dr - 


19 


i.e. 




= ^/X 2 


D = 


(56) 


The heavy quark diffusion coefficient D can be related to the momentum diffu¬ 
sion coefficient a and drag coefficient 77 through the Einstein relation 

2T2 _ T 
K Mrj 

It is used to describe the Brownian motion of a heavy quark in the hot medium and 
can also be related to the ratio of shear viscosity over entropy density 7^/5 liMI 
The shear ( 77 ) and bulk {Q viscosities of strong-interaction matter are extracted 
from correlation functions of the energy-momentum tensor, — 

I.e. using the current J^^(r, x) = T^^(t,x) in Eq. 




50 


one obtains 


Ti = IT lim hm - 

p^O 

C = — lim hm 

9 u;^0p^0 

k,l 


UJ 

pkk,ll(^^^ p) 


(57) 

(58) 


At present only a few calculations of transport coefficients have been performed 
in QCD with dynamical quark degrees of freedom. Most calculations utilized the 
quenched approximation at vanishing baryon number density.^I^ In the following 
subsections we will summarize current results for the transport coefficients intro¬ 
duced above. Eor further details on the calculation of transport coefficients on the 
lattice and earlier results see, for instance. Ref. 
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6.1. Electrical conductivity, charge diffusion and dilepton rates 


Probably the best analyzed spectral function is that of the light quark vector current 
which gives access to the electrical conductivity, cr(T) (Eq.[5^, soft photon emission 
rates (Eq. [5^ and thermal dilepton rates (Eq. [54| ) as well as the electric charge 
diffusion constant Dq = (jj^. The vector spectral function pv{^) = py(<^,0,T) 
has been calculated first in quenched QCEP^ using the MEM approaclUMIlll in 
which the analysis had only been constrained at large frequencies through a default 
model based on the free fermion spectral function and the largest lattices used only 
had = 16 points in the temporal direction. It has, however, been noticecP^ 
that this is not sufficient to get access to transport coefficients, which require the 
spectral functions at low frequencies to be linear in uj. Eurthermore, a larger time 
extenfP^ is needed to control the low frequency part of the spectral functions. 

Eirst continuum extrapolated results for the vector spectral function were ob¬ 
tained in quenched QCD using clover-improved Wilson fermionJ^^ at T 1.4 Tc. 
These results have been extended to three different temperatures recent 
In these calculations the thermal spectral functions are obtained by fitting the 
continuum-extrapolated two-point Euclidean correlation functions to an ansatz for 
the spectral function that is motivated by kinetic theory and perturbation theory 


^model 


(^) = Xq^BW 


ujT 






(59) 


Here Xq is similar to X 2 appearing in Eq. 55 but for light quarks and it is obtained 


by summing the zeroth component of the light temporal vector correlation function 
over space-time coordinates. The unknown parameters are thus the amplitude cbw 
and the width T of the Breit-Wigner function as well as the parameter k, which 
parametrizes deviations from the free theory behavior at large uj. The thermal 
spectral functions obtained in this way and the resulting thermal dilepton rates 
are shown in Eig. 19 In the high frequency region, i.e. for uj/T > 7, all the 


thermal spectral functions at the three temperatures above Tc are well described 
by leading order perturbation theory (Born rate) as well as Hard Thermal Loop 
(HTL) calculations. At small frequencies the lattice QCD results are significantly 
enhanced over the Born rate, indicating the presence of a transport peak, but are 
smaller than the HTL rates. In fact, the latter rises too rapidly at low frequency 
and would give rise to a divergent electrical conductivity. 


The finite intercept of p{uj)/{ujT) at u; = 0 seen in Eig. 19 gives the electrical 
conductivity (see Eq. 52), a/{Ce^T) G [0.21, 0.44].In these calculations a/T 
does not show any significant temperature dependency which may be due to the use 
of the quenched approximation. The results are consistent with a calculation that 
used staggered fermions, but no renormalized vector currents,^^^cr/ (Cgm^) — 0-4(1) 
for 1.5 < T/Tc ^ 2.25. These results are, however, substantially smaller than earlier 
calculations performed on coarse lattices.l^^ 


'This is also reflected in the temperature dependence of the temporal correlation functions divided 
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Fig. 19. (Left) Thermal spectral function at 1.1 Tc calculated in quenched QCD.EMl (Right) 
Temperature depe ndence of the quark-antiquark annihilation contribution to the thermal dilepton 
rate above Tc^^ 
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Fig. 20. (Left) Lattice QCD results for the e lectrical con ductivity cr/(CemT) as function of T/Tc 
including results obtained in quenched QCD^£E22E111 2-flavor QCD with ttItt — 270 MeW^ 
and (2+l)-flavor QCD with ttitt — 380MeV.E^ Note that the values of Tc obtained in these 
calculations are different, i.e. T^^ ° ~ 270 MeV, T^^ ^ ~ 208 MeV and T^^ ~ 185 MeV. 

(Right) Temperature dependence of charge diffusion coefficient 2ttTDq = 2'KTa j-y^. The results 
are obtained in (2+l)-ffavor QCD calculations on anisotropic lattices at one value of th e spa tial 
lattice cut-off, a, ~ 0.12 fm and with light quark masses corresponding to ttitt — 380 MeV.^^ 


Recently the calculation of the electrical conductivity has also been performed 
in QCD with dynamical quarks. The first calculation was done in 2-flavor QCD 
using 0{a) improved Wilson quarks corresponding to rrij^ ^ 270 MeV. The electrical 
conductivity of the quark-gluon plasma at T 250 MeV was found to be similar to 
the quenched result, cr/((7em^) = 0.40(12). 

First results from a calculation in (2-1-1)-flavor QCD performed on anisotropic 
lattices with quark masses corresponding to a pion mass 380 MeVEllE13 are 

shown in Fig. The simulation is performed at one value of the lattice cut-off. 

by which are temperature independent for this temperature range.EMl 
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Varying the temporal extent of the lattice, from between 48 and 16 a tempera¬ 
ture range from 0.63 Tc to 1.90 Tc is covered. The left hand panel of this figure shows 
temperature dependence of cr/((7emT) (square points). It stays constant within er¬ 
rors in the transition region and increases with temperature for T > T^. Also shown 
in the left hand panel of Fig. are electrical conductivities obtained in quenched 
QCEEMUnilllland 2-flavor QCD.S^The right hand panel of Fig. ^ shows the tem¬ 
perature dependence of charge diffusion coefficient Dq multiplied by 27rT. 27tTDq 
shows a dip near Tc which arises from the rapid exponential drop of the electric 
charge susceptibility at low temperature where the electrical conductivity of a pion 
gaJ^^ remains nonzero and vanishes only like \/T- 


6.2. Heavy quark diffusion 

Contrary to earlier expectations it has been found in heavy ion collision experiments 
at RHICll52ilI2l and LHCP^ that heavy-quark mesons show a substantial elliptic 
flow that is comparable to that of light-quark mesons. Moreover, also somewhat 
unexpectedly, heavy quarks are found to lose a significant amount of energy while 
traversing through the QGP. The latter has been called the heavy quark energy loss 
puzzle. Phenomenological explanations of these phenomena try to model the heavy 
quark diffusion in a hot and dense medium. This requires knowledge about the 
heavy quark diffusion coefficients 2ji54|i72|r7^ which can be determined in lattice 
QCD calculations. 

The charm-quark diffusion coefficient has been calculated in quenched QCD at 
three temperatures in the deconfined phas^^^ using lattices with temporal extent 
Nr = 48, 32 and 24 at a fixed value of the cut-off. Results from this calculation are 
shown in Fig.|^ The left hand panel shows the transport peak in the charmonium 
vector spectral function at three different temperatures above Tc. These transport 
peaks are obtained using a MEM analysis where an ansatz of Lorentzian form is 
used as the prior information for the low frequency part of the spectral function, i.e. 
ujr]l{uj^ The intercept at cj/T = 0 gives the charm quark diffusion coefficient 

(c.f. Eq. 55). It is found that 2itTD 2 and, with current understanding of 


uncertainties, this product is independent of temperature. These values of D are 
much smaller than the results from perturbative calculationi^^ and are very close 
to the value obtained from the AdS/CET correspondence.!^^ These different results 
are compared in the right hand panel of Eig. 

A direct determination of the diffusion coefficient for bottom quarks from bot- 
tomonium correlation functions is difficult as it would require quite small lattice 
spacings to accommodate the bottomonium states on the lattice. The method of 
choice here is to start with the infinite quark mass limit where one can use Heavy 
Quark Effective Theory (HQEET) to compute the diffusion coefficient of a static 
quark. In this approach the propagation of a heavy quark carrying color charge 
and its response to a colored Lorentz force can be described using linear response 
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Fig. 21. (Left) The low frequency transport part of the charmonium vector spectral func tion 
obtained from a MEM analysis of vector correlation functions in quenched lattice QCD.^^ The 
central lines are mean values while the bands reflect the statistical uncertainties. (Right) Temper¬ 
ature dependence of heavy quark diffusion coefficients multiplied by 2ttT in pure gluonic mat ter. 
The charm diffusion coefficients are obtained from the spectral functions shown in the left plot.^^ 
Also shown are results for t he stati c quark diffusion coefficient obtained by using the lattice dis¬ 
cretized versions of HQEFTEZEIlThe boxes show the statistical error while the error bars reffect 
the systematic uncertainties. The horizontal d otted line labels the value of 2ttTD in the heavy 
quark limit from the AdS/CFT correspondence and the short, horizontal soli d line indicates 
the value of 2ttTD from next-to-leading order pQCD calculations at cxg — 0.2.E^ 
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theory. This leads to the spectral analysis of a “color-electric correlator”[lIHE13 

^ l^{ReTv[U{l/T-,T)gEi{T,0)U{T-,0)gEi{0,0)]) 

= -(ReIY[C/(l/r;0)])-’ 

where U (r 2 ; ri) is a Wilson line in the Euclidean time direction and gEi^r^ x) denotes 
the color-electric field at time-space (r, x). The momentum diffusion coefficient e 
can then be extracted from the slope of the corresponding spectral function in the 
limit of vanishing frequency cj, 

^/t 3 = lim (61) 

Here Pe{^) is the corresponding spectral function that is related to Ge{t) via 
Eq. In the non-relativistic limit, i.e. for a heavy quark mass M ^ ttT, the 
momentum diffusion coefficient k is related to the (spatial) heavy quark diffusion 
coefficient D via Eq. 

The color-electric correlator Ge{t) is defined in terms of gluonic observables 
only. It thus suffers from a small signal to noise ratio, a problem similar to what 
one encounters in the calculation of viscositie^^^ (see next subsection). In both 
cases efficient noise reduction techniques are needed for a calculation of the gluonic 
correlation functions. In a recent stud}^^ the correlator Ge{t) has been calculated 
using multi-level update£SIE121 as well as the link integration technique EHEMl The 
improved correlator Gimpi^T) obtained this way at T 1.4 Tc for four different 
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Fig. 22. (Left) Continuum extrapolation of improved color-electric correlator G£;(r) = Gimp{'T) 
(see Eq. |60| normalized to the high temperature, free field theory result Gnorm{T) and the 
compariso n w ith NNLO and NLO perturbative QCD calculations at 1.4 Tc on quenched lat¬ 
tice QCD.[^3 (Right) fif fo fhe continuum extrapolated GimpIGnorm using an ansatz 
p'model(^^^ _ jcnajc {ApNNLo{(^) F Also shown are the separate contri buti ons to the 

correlator arising from the different parts of the fit ansatz for the spectral function.^3 


values of the lattice cut-off in quenched QCD has been extrapolated to the contin¬ 
uum limit. Results are shown in the left hand panel of Fig. In order to extract 
the momentum diffusion coefficient from this correlator an ansatz for the spectral 
function has been used.^^^ The corresponding fit is shown in the right hand panel 
of Fig. It yields, 

i^/T^ = 2.5(4) at T - 1.4 T^. (62) 

This corresponds to 27tTD = 5.0(8) which is consistent with similar lattice QCD 
studies using HQEFT.Id^ However, it is a factor 2 to 3 larger than the charm quark 
diffusion coefficient discussed abovd^^ (see Figj^. 

The heavy quark diffusion coefficient D has also been calculated in perturbative 
QCD in both leading and next-to-leading order as well as from the AdS/CFT corre¬ 
spondence. At as ~ 0.2 the leading order pQCD calculatiorP^ gives 27tTD 71.2 
while 27tTD ^ 8.4 is obtained in next-to-leading order calculationsIn the strong 
coupling limit 27tTD = 1 is obtained from the AdS/CFT correspondence^ (see 
Fig^. 


6.3. Shear & bulk viscosities 


The computation of shear and bulk viscosities has been attempted already in the 
80 ’sJIHni However, it turned out that the signal to noise ratio in calculating Euclidean 
energy-momentum correlation functions is very small and high statistics or improved 
algorithms are needed. The latter became available for quenched QCD calculations 
with the multi-level algorithm^^^ which have been used for the calculation of shear 
and bulk viscosities EHESZl Results from such an analysis by Meyei^i^ and those 
with direct high statistics calculations by Nakamura & SakaP^are shown in Eig. 23 
















THERMODYNAMICS OF STRONG-INTERACTION MATTER FROM LATTICE QCD 43 


The ratio obtained in the former calculation are smaller than those obtained in the 
earlier direct simulations; rj/s = 0.102(56) at 1.24 Tc and 0.134(33) at 1.65 Tc 
obtained in quenched QCD on lattices with temporal extent = 8 .^^^ This shows 
that noise reduction techniques like the multi-level algorithm are mandatory for a 
successful calculation of viscosities. Results on the bulk viscosity to entropy ratio 
C/s show that this ratio rapidly becomes small above At T > 1.2 Tc it is 

smaller than 77 /s and, in fact, within errors it is consistent with zeroEHEMllHl 
Still the calculations of viscosities are performed on lattices with rather small 
temporal extent compared to those used in calculations of the electrical conductivity 
and diffusion constants. Systematic uncertainties in these calculations need to be 
better controlled in future. Here it will also be helpful to combine lattice QCD 
calculations with information from analytic approaches that put constraints on the 
structure of spectral functions, e.g. QCD sum rules 
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Fig. 23. (Left) Temperature dependence of th e shear viscosity, 77/5, of gluonic matter obtaine d by 
using noise-reduction techniques (squares^P^ and direct high statistical calculations (circlesp^ 
on Nr = 8 lattices. The dott ed line corresponds to the conjectured lower bound of 77 /s from 
AdS/CFT correspondence(Right) Temperature dependence of the bulk viscosity, (^/s^ of 
gluonic matter.^^ The square points include the statistical uncertainties while the solid black 
bars denote the systematic uncertainties. Results for the same quantity from Ref. |188| which are 
not shown are consistent with zero at T G [1.4 Tc, 1.8 Tc]. Data are taken from Refs. |181|188|18^ 


6.4. Transport coefficients of second order hydrodynamics 

Due to the importance of viscous effects in the evolution of hot and dense matter 
created in heavy ion collision there is considerable interest in extending the hydrody¬ 
namic modeling beyond leading order gradient expansions of the energy momentum 
tensor. A second order gradient expansion is parameterized by additional transport 
coefficients,^!^ which may become accessible to lattice QCD calculations liMriHH 
One of these new transport coefficient^ is which controls the momentum depen- 


^Although this particular second order coefficient is also called k it should be noted that it is not 
related to the momentum diffusion coefficient k introduced in Eq. 
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dence of space-like correlations of the energy momentum tensor. It is defined as 
the leading order coefficient of the Taylor expansion of the retarded correlator of 
energy-momentum tensor in momentum p about zeroflMIlZI 

/ oo 

dt dx e*P’'([Ti 2 (T,x),T i 2(0, O)])0(i) 

-OO 

= G^(0) + f |p|^ + 0(|ph. (63) 

Since the difference between the retarded correlator and the Euclidean correlator 
at vanishing frequency cj = 0 is just a contact term, which has been shown to be 
momentum independent k can directly be extracted from the corresponding 
Euclidean correlator, 


K = 


2 lim 

IpKo 


dG^{uj = 0, p) 

d|p|2 


(64) 


This is considerably simpler to compute than first order coefficients such as the shear 
or bulk viscosity as it does not require analytic continuation from imaginary time 
to real time. The only complication is that large spatial lattice sizes are required 
to get access to small momenta in the lattice QCD calculation. 



Fig. 24. Transport coefficient of second order hy drody namics n/T‘^ as a function of T. The l ines 
denote the results from AdS/CFT correspondencand lattice perturbation theory (LPT),E^ 
respectively. 


The second order hydro-coefficient k has been computed in quenched QCD on 
lattices with temporal extent Nj- = 8 and 6 and a large spatial size.^^^ Results from 
this calculation are shown in Eig. 24 It yields k/T‘^ = 0.36(15). Within the current 
large errors no further temperature dependence is observed in the temperature range 
2Tc <T < lOTc and results are compatible with both perturbative QCEEMHHand 
AdS/CET resnlts.l^ 
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7. Open heavy flavors and heavy quarkonia 


Owing to their large masses heavy flavors, such as charm and bottom, are produced 
only during the earliest stages of a heavy ion collision and get affected by the hot and 
dense medium during its entire evolution. They then emerge as hadrons with open 
heavy flavors or as quarkonia. Thus, heavy flavors provide unique identifiable probes 
against which the temperature and coupling of the medium created during heavy- 
ion collisions can be calibrated. The melting of charmonium states has long been 
considered as an important signature for the formation of quark-gluon plasma.^^ 
Understanding its production rates and abundance can provide detailed information 
on properties of the matter formed in a heavy ion collision. We will discuss current 
results on the temperature dependence of quarkonium bound states in subsection 


7.2 In addition, understanding the fate of open flavor bound states also is of 
importance for the understanding of the strongly interacting regime in the QGP 
close to the transition temperature. The question whether or not heavy flavor 
mesons can survive above plays a crucial role in the phenomenological modeling 
of the heavy quark energy loss.l^^®^^ Lattice QCD results on the dissociation of 
open charm hadrons will be discussed in the next subsection. 


7.1. Melting and abundances of open charm hadrons 


Similar to the case of strangeness,^ discussed in Section 3.2, the deconfinement or 


melting of open charm hadrons can be probecP^ by using fluctuations of the charm 
quantum number (C) and its correlations with baryon number (B) fluctuations. 
To this end one may consider appropriately chosen ratios of correlations between 
net-baryon and net-charm number fluctuations. One such ratio is, for instance, 
X2‘f lx Since at low temperatures, for a gas of uncorrelated charmed hadrons, 

the thermodynamic contributions of the much heavier \C\ = 2,3 baryons can be 
safely neglected, this particular ratio is dominated by \C\ = 1 baryons and thus 
yields X 2 ‘f IXi^ = 1^1 = 1- fhe other hand, at high temperatures, for a free 
gas of massive charm quarks with masses much larger than T, the above ratio 
= 1^1 = 1/3- Thus, this quantity can be viewed as a measure for 
the baryon number associated with the predominant degrees of freedom that are 
carriers of the charm quantum number in a given temperature region. Its behavior 
can indicate the liberation of quark-like degrees of freedom with fractional baryon 
number. 

Fig. 25 (left]P^ shows that close to the chiral crossover temperature, Tc = 
154(9) MeV, the quantity X 2 ‘f IXi^ starts deviating from unity. This indicates 
that beyond this temperature fractionally charged quark-like degrees of freedom 
start appearing and open charm hadrons start to melt. The left hand panel of 


SRere one makes use of the fact that the Boltzmann gas approximation is still a good approximation 
for a charm quark gas at a few times the transition temperature. This is the temperature range 
of interest to heavy ion collision experiments. 
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Fig. also shows a comparison with similar quantities involving both baryon- 
strangeness correlations as well as baryon-electric charge correlations. The first of 
these quantities is sensitive to the deconfinement of strange quarks and the second 
one is also sensitive to the deconfinement of light up and down quarks. It is clear 
from the figure that the onset of deconfinement for the charm, strange as well as 
the up/down quarks happens in the same chiral crossover region. 
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Fig. 25. (Left) Ratios between cumulants of correl ation s of net baryon number with net charm, 
net strangeness and net electric charge fluctuationsDeviations of these ratios from unity in 
the crossover region (shaded band) suggest that the onset of the melting/deconflnement of open 
charm hadrons as well as open strange hadrons and hadrons with light up/down quarks starts in 
the same temperature range (see text for details). (Right) Thermodynamic contributions of all 
charmed baryons (top), all charged charmed baryons (mi ddle) and all strange charmed baryons 
(bottom) relative to that of corresponding charmed mesonsThe dashed lines (PDG-HRG) are 
predictions for an uncorrelated hadron gas using only the PDG stateJ^ . The solid line s (QM- 
HRG) are similar HRG predictions including also the states predicted by the quark modelIl2lIl2ll 
The dotted lines (QM-HRG-3) are the same QM predictions, but only including states having 
masses less than 3 GeV. 


The right hand panel of Fig. shows ratios of generalized susceptibilities that 
are constructed such that in a gas of uncorrelated hadrons the numerator would 
correspond to the partial pressure of baryons with quantum numbers selected by 
the other indices, while the denominator would correspond to the partial pressure of 
corresponding mesons. E.g. Xi^ will give the partial pressure of charmed baryons 
and Xa ~ Xis^ filters out the corresponding partial pressure of charmed mesons. 
The ratios x?f/{xi - xf,'"). xf,T/(x?F - xf^T) -tid - xfi^) 

thus correspond to ratios of partial pressures arising from charmed baryons to that 
from open charm mesons, charged-charmed baryons to open charm charged mesons 
and strange-charmed baryons to strange-charmed mesons, respectively. As can be 
seen in Fig. (right) results for these ratios at low temperatures are well described 
by a hadron gas that uses open charm resonances obtained in quark model calcu¬ 
lations (QM-HRG) while, not surprisingly, it differs strongly from a simple hadron 
resonance gas based only on the very few experimentally known charmed hadrons 
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(PDG-HRG). This observation is similar to that made for strange hadron fluctu¬ 
ations in Section]^ (right panel of Fig. 18). These observations clearly provide 


evidence for contributions from experimentally yet unobserved charmed/strange 
hadrons to charm and strangeness fluctuations in the vicinity of the crossover tran¬ 
sition temperature. In additions, the agreement between lattice QCD data and QM- 
HRG model calculations below Tc and the onset of deviations just in the transition 
region also suggests that open charm hadrons start to dissociate in the crossover 
region. This conclusion is consistent with that drawn from the temperature depen¬ 
dence of the ratio in the left hand panel of Fig. which has the advantage of being 
independent of details of the charmed hadron spectrum. 


7.2. Heavy Quarkonia 

Inside the QGP the interaction between a heavy quark anti-quark pair gets weak¬ 
ened due to the screening effects of the intervening deconfined colored medium. 
This telltale signature of the presence of a color deconfined medium is expected to 
be manifested through melting of heavy quarkonium states,!^ i.e. bound states of 
a quark anti-quark pair. Melting of heavy quarkonium states in the QGP is also 
expected to follow a distinctive sequential pattern with the smallest, most tightly 
bound quarkonium state surviving up to the highest temperature, effectively serv¬ 
ing as a ‘thermometer’ to probe the temperature of the medium created at RHIC 
(for a recent review see Ref. 206). Since quarkonia do not carry flavor quantum 


numbers, unlike the open heavy flavor hadrons, their melting cannot be accessed 
in lattice QGD calculations through quantum number correlation based studies as 
discussed in Section |7.I[ Gurrently, at least three approaches are being actively 
pursued to study thermal modifications of heavy quarkonia on the lattice. The 
first approach is to use potential models with heavy quark potentials computed on 
the lattice. The second approach relies on the extraction of quarkonium spectral 
functions from Euclidean temporal correlators. The third approach involves spa¬ 
tial correlation functions of quarkonia and the study of their in-medium screening 
properties. Here, we summarize the recent progresses and developments concerning 
these three approaches. More details about earlier lattice results can also be found 
in Refs. [206lMi8l 


7.2.1. Heavy quark potential 

Properties of heavy quarkonia in the vacuum have been successfully described by the 
potential model approach. In this approach, the interaction between a heavy quark 
pair forming the quarkonium is described by an instantaneous potentialE2lI2l3 Due 
to its success at zero temperature, the potential model approach has also been ap¬ 
plied at nonzero temperature, by making the potential between the heavy quarks 
temperature dependent The temperature dependent potentials used in these 
calculations are based either on model calculations or on finite temperature lattice 
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QCD results for the free energy and the internal energy. However, justification 
for the use of these potentials in a Schrodinger equation is mostly based on phe¬ 
nomenological arguments and has no rigorous connections to the real-time evolution 
of heavy quarkonia described by the Schrodinger equation. 

Recently, a non-relativistic effective theory approach at nonzero temperature, 
relying on a particular hierarchy of the relevant scales, has been developed 1^2^214] 
Generally by integrating out the hard energy scale, i.e. the rest mass of the heavy 
quark m, Non-Relativistic QCD (NRQCD) effective theor}!!!!!^!^! is obtained, and 
by further integrating out the soft scale, i.e. the typical momentum exchange be¬ 
tween the bounded quarks of 0(mn), so-called potential non-relativistic QCD (pN- 
RQCD), is obtained.In pNRQCD a heavy quark bound state can be described 
by a two-point function satisfying a Schrodinger equation 

idtD>^{t,f^ = (65) 


where D^^(t^r) is the real-time forward heavy quark pair correlation function in 
the non-relativistic limit. The potential V(t^r) in such a Schrodinger equation is 
well defined and turns out to be complex valuedIts real part reflects color 
Debye screening effects while the imaginary part is related to Landau-damping, i.e. 
the scattering of quarks with the constituents of the medium and the absorption 
of gluons from the medium via singlet-octet transit ion.It has been shown in 
Ref. 218 that in the limit m ^ oc the leading part of the potential V{t^r) can be 
obtained as 


V(r) = lim 

t^oo 


idtW{t,r) 

W{t,r) 


( 66 ) 


where IT(t, r) is the real-time thermal Wilson loop. Through a Fourier transforma¬ 
tion one obtains. 


/ + 00 

dcj r) 

-oo 




(67) 


where p{uj,r) is the spectral function of the real-time Wilson loop, 
equation can be analytically continued to the Euclidean time, giving 

^+oo 

duj p{uj^ r) . 


r 

W(r,r)= / 

J —( 


The above 


( 68 ) 


In Eq. 68 the Wilson loop, W (r, r) is the usual Euclidean-time Wilson loop which 
can be obtained from lattice QCD calculations. As if there exists a well defined 
lowest lying peak structure in the spectral function p{ijj^ r) it would dominate the 
dynamics in the Wilson loop in the late time limit thus giving to the complete 
information of the heavy quark potential as shown in Eq.|^ However, as argued in 
Ref. |219| the late time dynamics is not completely separated from the early time non¬ 
potential physics and to take these effects into account a skewed Lorentzian function 
including the early time contribution for p(co’,r) is required to fit the Wilson loop. 
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Fig. 26. Real (left) and imaginary (right) parts of the heavy quark potential calculated in 
quenched lattice QCD.I^^ Values are shifted for better visibility. Gray circles denote color singlet 
free energies and solid lines are leading order hard thermal loop results. 


The position and width of this skewed Lorentzian peak are connected to the real and 
the imaginary parts of the potential, respectively. Thus, by extracting the spectral 
function, p(co’,r), from the Euclidean-time Wilson loop|^ lT(r, r), obtained from 
lattice QCD calculations and by subsequently fitting this p(co’,r) with a skewed 
Lorentzian ansatz the real and the imaginary part of the heavy quark potential, 
E(r), can be obtained. 

This procedure has been tested using the Wilson loop calculated using hard 


shows the real (left) and imaginary (right) parts of the heavy quark potential in the 
pure gluonic medium .1^^ The real part of the heavy quark potential turns out to be 
close to the heavy quark singlet free energy, defined through the spatial correlation 
function of the Polyakov loop and its conjugate,^^^ at all temperatures. At high 
temperatures the imaginary part of the heavy quark potential is close to the leading 
order hard thermal loop result, but at low temperatures it lies below hard thermal 
loop result. 

Using the same procedure the real part of the heavy quark potential has also been 
extracted in (2-f-l)-flavor QCEP^on lattices with temporal extent of only = 12. 
Of course, on such small lattices it is hardly possible to extract reliable information 
on the imaginary part of the heavy quark potential. Nonetheless, similar to the 
quenched limit result the real part of the heavy quark potential is found to be close to 
the singlet free energies. However, string breaking is not observed at T < Tc ~ 174 
MeV up to ^ 1.2 fm. This could be due to the large pion mass (tUt^ ~ 300 MeV) 
employed in this exploratory QCD calculation. 

The heavy quark potential has also been extracted from the Wilson line correla¬ 
tion function calculated on 48^ x 12 lattices at T = 250 MeV and 305 MeV using a 
fit ansatz motivated by hard thermal loop calculations.!^^ These calculations used 

^As illustrated in Ref. |22Q| in order to avoid the cusp divergence and the large noise to signal ratio 
in computations of Wilson loops, some other observables, e.g. Wilson line correlation functions 
defined in the Coulomb gauge, can be a substitution and are actually used in lattice calculations. 


thermal loop perturbation theor 3 @^as well as for the pure gauge theory.!^^ Fig. 26 
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HISQ action with quark masses that would correspond to ~ 160 MeV in the 
continuum limit. In this calculation the real part of the potential was found to 
be equal or larger than the singlet free energies at these two temperatures and it 
always was smaller than the zero temperature potential. Moreover, the imaginary 
part is found to be of similar size as predicted by leading order hard thermal loop 
perturbation theory.l^^ 


7.2.2. Spectral functions of quarkonia 


The spectral function of quarkonium can be extracted from lattice QCD calculations 
of the two-point correlation functions of the heavy quark currents in the Euclidean 
time 

G{t,T) = J (fx (q{T, x) r g(r, x)^ (^^(0,0) T g(0,0)) ^ ^ , (69) 

where q represents the heavy quark field. T are Dirac matrices defining the spin 
structure of quarkonium, in particular T = 1 , 75 , 7 ^, 757 ^ corresponds to the scalar, 
pseudo-scalar, vector and axial-vector quarkonium states, respectively. Signatures 
of medium modification and melting of quarkonium are reflected in the structure of 
the spectral density p{uj^T). It can be obtained from the quarkonium correlator in 
the Euclidean-time 



n( ■ T'l ~ 

27r^^ ’ ^ sinh(w/2T) 


(70) 


While p{uj,T) is a continuous function of the frequency uj over an infinite range, 
G{r,T) is calculable on the lattice only at a finite number of discrete points along 
the temporal direction due to the finiteness of the lattice. As has been discussed in 
Section [^inverting Eq. 70 to get spectral function p(co’, T) from the correlation func¬ 
tion G(r, T) is a typical ill-posed problem and a suitable inversion method is needed 
to solve the problem. The commonly used inversion methods for the extractions of 
quarkonium spectral functions is the Maximum Entropy Method (MEM).Vari¬ 
ants of this have recently been suggested which include a modified MEM with an 
extended search spacdHHilHH^ and a Bayesian method which is analogous to MEM 
but uses a different prior distribution.!^^ 

The charmonium spectral functions at finite temperature have been studied 
extensively for quenched lattice QCD.Recently quenched lattice QCD calcu- 
lation^^^ have been performed using lattices with large temporal extents, leading 
to quite reliable extractions of the charmonium spectral functions. Eig. illus¬ 
trates the spectral functions of J/T (left) and pc (right) for various temperatures. 
These studies suggest that all charmonium states are dissociated for T > 1.5Tc in a 


gluonic plasma.l^^ Investigations on charmonium spectral functions in lattice QCD 
calculations with dynamical quarks lead to consistent conclusions EMESI 

Lattice QCD calculations for bottom quarks are technically challenging since 
present computational limitations do not permit calculations with sufficiently fine 
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Fig. 27. Temperature dependence of J/T (left) and r]c (right) spectral functions obtained from 
quenched lattice QCD calculations.!^^ 


lattice spacings that are much smaller than the inverse of the heavy bottom quark 
mass. Large discretization errors arise from the large heavy-quark bare mass in 
units of the lattice spacing, m^a. Nonetheless, although it is computationally very 
demanding relativistic bottomonium correlations functions have been computed re¬ 
cently in quenched QCD using isotropic lattices with A^cr=192 and large temporal 
extent,!^^ = 96 and 48 at the same value of the lattice spacing, a 0.01 fm, 
i.e. miya 0.2. The temporal lattice sizes correspond to temperatures T = 0.7 Tc 
and 1.4 Tc. In this calculation it has been found that at 1.4 Tc only the bottomo¬ 
nium correlation function in the vector channel show significantly smaller thermal 
modifications than the charmonium correlation function in the same channel.l^^ 
To circumvent the problem of large cut-off effects in calculations with bottom 
quarks without performing calculations at extremely small lattice spacings, one 
needs to adapt a lattice discretization scheme for heavy quarks that is capable 
to describe heavy quark physics in some chosen kinematic regime also at mod¬ 
erately small values of the cut-off. One such approach is Non-Relativistic QCD 
(NRQCDjSmSlll where the heavy quark mass is assumed to be much larger than 
the inverse lattice spacing, but the momentum dependence of the heavy quark en¬ 
ergy is included in the non-relativistic limit. NRQCD does not possess a proper 
continuum limit because the radiative corrections to coefficients of the NRQCD La- 
grangian diverge as m^a ^ 0. However, since the non-relativistic quark field is not 
compactified along the Euclidean-time direction, the relation between Euclidean¬ 
time quarkonium correlation functions and the corresponding spectral functions 
becomes simpler, 

G(T,T) = ^°°^p(a;,T)e-“V (71) 

giving rise to a temperature independent integration kernel, exp(—c^r). Moreover, 
the transport peak is also absent in the NRQCD spectral function, which makes its 
determination much easier. 
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uj (GeV) 



uj (GeV) 


Fig. 28. Temperature dependence of spectral functions for the T (top) and Xhi states (bottom). 
Results are obtained by using a lattice discretized form of NRQCD.l^^ 


Treating the bottom quark within the NRQCD framework the fate of bottomo- 
nium states has been studied on anisotropic lattices using Nf = 2 and 2 + 1 of 
improved Wilson quarks with ~ 400 MeV and ~ 500 MeV, respec- 

tivelyUmmi Fig. 


28 


shows the temperature dependence of the T (top) and 
Xhi (bottom) spectral functions extracted using the This indicates that 

the S-wave ground state T survives up to 1.9 Tc while the P-wave ground state 
Xhi melts just above However, a recent lattice QCD stud}P^ that uses a 

different Bayesian method suggested in Ref. 226 for the extraction of the spectral 
functions finds that the P-wave ground state Xhi i^ay survive up to 1.6 Tc. This 
analysis, however, has been performed on lattices with a rather small temporal 
extents, W = 12. 
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7.2.3. Spatial correlation functions of quarkonia 


Similarly as the temporal correlation functions spatial quarkonium correlation func¬ 
tions are defined as the sum of heavy quark current-current correlation functions 
instead over a ‘funny’ space of r, x and y 

Giz,T) = J dr J dxdy (^(q{T,x)rq{T,x)j(^q{0,6)'['q{0,6)y^ , (72) 

are also related to the quarkonium spectral functions as follows 


G{z,T) = 


7-00 27r 7o 


dct; p{u:,pz,T) 
27r UJ 


(73) 


The spatial correlation functions, however, are sensitive also to the momentum de¬ 
pendence of spectral functions which is reflected in an additional integration over 
the momentum pz along the 2 )-direction. Despite of its more complicated connec¬ 
tion to the spectral function compared to that of temporal correlators (see Eq. 
70), the spatial correlation function has some distinct advantages. Since the spa¬ 


tial separation is not limited by the inverse temperature, the spatial correlation 
function can be studied at separations larger than 1/T and thus is more sensitive 
to in-medium modifications of quarkonia. Furthermore, in contrast to the tem¬ 
poral correlation functions, the spatial correlation functions can be directly com¬ 
pared to the corresponding vacuum correlation function to quantify modifications 
of the in-medium spectral function as the entire temperature dependence of spa¬ 
tial correlation functions emerges from the temperature dependence of the spectral 
function (Eq. 73), while in the case of temporal correlation functions (Eq. 70) this 
temperature dependence is folded with the temperature dependence of the kernel, 
cosh {uj (r — 1/2T)) / sinh {uj/2T). 

While the relation between spectral functions and spatial meson correlators is 
more involved, in some limiting cases it becomes more intuitive. At large distances 
the spatial correlation functions decay exponentially, G{z,T) ^ exp(—M(T) 2 :), 
characterized by the screening mass M. At low temperatures this exponential drop 
of a mesonic correlator projecting on a given quantum number channel E is con¬ 
trolled by the bound state with smallest mass in this channel. The spectral function 
then is dominated by this state. 


p{oj,pz, T) ~ <5(^2 -pI- ml), 


(74) 


where M becomes equal to the pole mass mo of the meson. At high temperatures, 
when the mesonic excitations have completely dissolved and become uncorrelated, 
the spatial meson correlation functions describe the propagation of a free quark 
anti-quark pair and the screening mass is given bjEH 

Mfree = fml + (7rT)2 + + (nT)^, (75) 

where m^^ and rUg^ are the masses of the quark and anti-quark that form the meson. 
This form of the screening mass in the non-interacting limit is a consequence of the 
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Fig. 29. (Top) Temperature dependence of ratios of spat ial correlation functions to those at zero 
temperature for D* (top left) and J/4^ (top right).^^^ (Bottom) Temperature dependence of 
screening masses of ground state sc mesons (left) (l^)? (0+), D* (1“) and Dg (0“) 

and ground state charmonia (right) Xci (l"^"^)? Xcq (O"^"^)? T/4^ (1 ) and r]c (0 I- ) 1245] q^i^g 

horizontal solid lines denote the corresponding zero-temperature meson masses and the dashed 
line depicts the non-interacting theory limit of a freely propagating quark anti-quark pair. The 
yellow band denotes the chiral crossover temperature region, i.e. Tc = 154 ± 9 MeV. 


fermionic (anti-periodic) boundary condition along the Euclidean-time direction, 
which leads to the appearance of a smallest non-zero Matsubara frequency, ttT, in 
the quark and anti-quark propagators. As the bosonic meson state is dissolved in 
the non-interacting limit the screening mass reflects two independently propagating 
fermionic degrees of freedom. Thus, the transition between these two limiting values 
for the screening mass can be used as an indicator for the thermal modification and 
eventual dissociation of charmonia and open charm mesons. 

A first analysis of spatial correlation functions of charmonia was performed 
using staggered fermions with the p4 action.^^^ Recently, a detailed comparative 
study of the spatial correlation functions of charmonia and open charm mesons was 
performed using the HISQ action.^^^ Fig. shows ratios of the spatial correlation 
functions for the (top left) and the J/T (top right) at different temperatures 
to the corresponding zero temperature results. As discussed before, such a ratio 
can directly probe the thermal modifications of the spectral functions. If there is 
no change in the meson spectral functions, these ratios will be equal to unity and 
deviations from unity will indicate in-medium modification of the meson spectral 
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functions at non-zero temperature. As can be seen, even for T ^ Tc ~ 154 MeV 
the ratio for the open charm sector (Fig. (top left)) shows significant deviation 
from unity. This is completely consistent with the previously discussed quantum 
number correlation based studies in Section [7T] which concluded that open charm 
mesons start to melt in the vicinity of On the other hand, for the charmonia 

the ratio of non-zero to zero temperature correlation functions remains very close 
to unity for T < 170 MeV and significant deviations, comparable to the open 
charm sector at T ^ Tc, only show up for T > 200 MeV (Fig. (top right)). 
This indicates that up to T 1.1 Tc J/T exists almost like a vacuum state and 
significant thermal modifications appears only for T > 1.3Tc. The same physics 
is also reflected in temperature dependence of the screening masses of open charm 
mesons and charmonia as shown in Fig. (bottom). The screening masses of all 
open charm mesons (bottom left) start to deviate from unity significantly already 
in the chiral crossover region. The J/T and r]c screening masses (bottom right), 
however, remain to be almost identical to their vacuum pole masses up to T 
1.1 Tc. The screening masses for Xco Xci deviate from their vacuum masses 
even below Tc indicating significant thermal modifications of these states. At higher 
temperatures, i.e. T > 2 Tc the screening masses in all channels increase linearly 
with temperature. This rise agrees well with that of a freely propagating quark 
anti-quark pair. 

8. QCD in external magnetic fields 

In recent past it has been realized that large magnetic fields created during the early 
stages of heavy-ion collisions may give rise to fascinating observable effects induced 
through the coupling between the magnetic field and the chiral anomalyESSIl This 
observation also motivated a plethora of activities including lattice QCD studies of 
hot-dense strong-interaction matter under the influence of external magnetic fields. 
This section presents a very brief summary of these lattice QCD studies. More 
comprehensive reviews on this topic can be found in Refs. |247f[25Q| 

An external magnetic field along the ^-direction, B = Bz, can be induced 
by choosing electro-magnetic gauge fields: = Bx and = 0. 

On the lattice this can be implemented simply by multiplying the SU{3) gauge 
field variables Un,ji with the corresponding U{1) phase factors: Un,y = qDux 
and Un,x = Lin^z = '^n,f = 1- Here, q is the (electric) charge of a quark 
and n = {nx^fiy^nz^rir) denotes a lattice site, with nx,y,z = 1 ••• and 
Ur = 1 ... Nx. This choice ensures that the magnetic fiux, a^T, is constant through 
all the plaquettes in the x-y plane, except at the boundary (N^j^ny^nz^rir) owing 
to the periodic boundary condition on the gauge fields along the spatial directions. 
Thus, to preserve the smoothness of the external magnetic field across the bound¬ 
ary and the gauge invariance of the fermion action the U{1) factor on the boundary 
links must be modified: 'n(Ar^,n qN^Buy ^ gy making such a choice 

one pays the price that the magnetic fiux becomes quantized, o?B = 27rnB/{qN^) 
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Fig. 30. (Left) The QCD crossover temperature as a function of the magnetic field obtained from 
the inflection point of the quark condensate (red band) and strange quark susceptibility (blue 
band).I^^ (Right) The magnetic susceptibility xb^ be. the correction t o th e pressure proportional 
to quadratic power of the magnetic field, as a function of temperature .1^^ 


with = 0 ... N‘^. Thus, at a given temperature T and on a finite volume there 
is a minimal value for the constant, external magnetic field B = 27rT‘^/{qN^). 
On the other hand, for given temporal extent, of the lattice the maximal value 
of the magnetic field achieved in this case is B = 27 tT‘^N^/ q. Most importantly, 
since the external magnetic field is induced by multiplying the SU{3) gauge field 
variables with purely imaginary phase factors the fermion determinant remains real 
and positive, allowing direct Monte-Carlo sampling. 

The QCD crossover transition in an external magnetic field has been extensively 
studied using lattice QCD. In contrast to the earlier effective model based predic¬ 
tions and the results of first lattice QCD calculationJSHEISl using unimproved stag¬ 
gered fermions with heavy pion masses, the QCD crossover temperature Tc turns 
out to decease with increasing magnetic field.^^^^M^^ Fig. 30 (left) illustrates the 
decrease of Tc as function of the magnetic field.The reduction of Tc follows from 
the reduction of the chiral condensate with increasing 5, a phenomenon known 
as the inverse magnetic catalysis, for the relevant intermediate temperature range. 
However, for low as well as high enough temperatures the chiral condensate shows 
the expected magnetic catalysis, i.e. an increase with increasing 

The response of the hot strong-interaction matter to an external magnetic field 
has also been investigated through lattice QCD calculations of the magnetic suscep¬ 
tibility. The magnetic susceptibility is the second derivative of the QCD partition 
function with respect to the external magnetic field B and, thus, gives the order B^ 
corrections to the pressure. Since the standard approach leads to discrete, quantized 
values for B on the lattice, unlike for the susceptibilities of the chemical potentials, 
taking direct derivatives of the QCD partition function with respect to B is not 
straightforward. Lattice calculations of the magnetic susceptibility have been per¬ 
formed using three different approaches: (i) by computing the free energy difference 
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through an integration over the varying magnetic field(ii) through com¬ 
putations of the magnetization via pressure anisotropies parallel and orthogonal to 
the direction of the magnetic field|H£llH^ (iii) by adopting a non-uniform external 
field such that the net magnetic fiux across the x-y plane of the whole lattice is zero, 
resulting in a non-quantized magnetic field and allowing for a direct computation of 


(right), all these different approaches produce more or less consistent results. The 
QGP phase turns out to be paramagnetic. However, for low enough temperatures 
QCD is weakly diamagnetic, in accordance with the hadron resonance gas model.l^^ 
These results also suggests that for T < 300 MeV the order magnetic contribu¬ 
tion to the QCD pressure is at best a few percent for a magnetic field of the order 
of 10^^ T, which is relevant for heavy-ion collision experiments.!^ 

Lattice QCD studies have also been carried out to investigate the chiral magnetic 
effect,!^ namely the phenomenon of (electric) charge separation along the direction 
of the magnetic field in presence of chiral anomaly induced topological fiuctuations. 
Various different avenues have been pursued by studying: (i) charge ffuctuations on 
given topological background(ii) enhanced ffuctuations of the electric current 
along the magnetic field,(iii) dependence of the electric current on a given axial 
chemical potential(iv) correlations of the electric polarization and the topolog¬ 
ical charge.!^^ In all cases, the chiral magnetic effect has been observed, however, 
in a suppressed magnitude compared to the model based expectations. Obviously, 
in lattice QCD studies of the chiral magnetic effect, which is closely related to the 
chiral anomaly of QCD, it would be of great interest to have lattice QCD calcu¬ 
lations perfomed with a chiral fermion formulation. Fully dynamical calculations 
with DWF or overlap fermions at present do not yet exist. 


the derivatives of the partition function with respect to As shown in Fig. 


9. Summary 

In this review we have discussed selected topics in finite temperature QCD that 
have a direct link to the ongoing experimental study of the phase diagram of strong- 
interaction matter and the properties of matter formed in heavy ion collisions. We 
naturally focused on topics to which lattice QCD calculations can contribute and 
have lead or may lead to definite, quantitative answers in the future. 

Lattice QCD calculations at finite temperature started more than 30 years ago 
with simulations of pure SU{N) gauge theories on computers that delivered less 
than one Mega-Flop/s peak-performance. Nowadays these calculations are per¬ 
formed with dynamical quarks and the correct mass spectrum of QCD on comput¬ 
ers that reach a peak-performance of more than ten Peta-Flop/s. The field thus 
could utilize computing resources with a peak-performance that increased by more 
than ten orders of magnitude, or doubled the speed almost every year. Although we 
know quite well from the physics of the hadron gas that an exponential rise cannot 
continue for ever, we still do not seem to have reached the critical point and we may 
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also in the future expect to gain further inside into the physics of strong-interaction 
matter through numerical lattice QCD calculations. 

During the last years lattice QCD calculations have accomplished two long¬ 
standing goals: (i) the pseudo-critical temperature of strong-interaction matter 
with physical light and strange quark masses has been determined at vanishing 
net-baryon number, Tc = (154 zb 9) MeV, and (ii) the equation of state has been 
calculated in a wide temperature range. It now is frequently used in the hydrody¬ 
namic modeling of heavy ion collisions. The extension of these results to non-zero 
chemical potential is well controlled at least in leading order Taylor expansion of 
the pressure in terms of baryon, electric charge and strangeness chemical poten¬ 
tials and we will soon have results for bulk thermodynamics in a range of baryon 
chemical potentials /iB/T<3, which will cover most of the parameter space that 
can be explored in heavy ion experiments at RHIC and the LHC. These Taylor ex¬ 
pansions also provide specific information on fluctuations of and correlations among 
conserved charges. Lattice QCD calculations of these observables have reached a 
stage where guidance can be provided to the experimental search for the critical end 
point and, with some care, a comparison between experimental results for measured 
proton number, strangeness and, in particular, electric charge fluctuations on the 
one hand and lattice QCD calculations on the other hand starts to become possible. 

A major focus in the ongoing experimental heavy ion program is to explore in 
detail the properties of matter formed in heavy ion collisions. The analysis of trans¬ 
port coefficients, thermal masses and screening lengths does play an important role 
in this effort. Lattice QCD calculations of these quantities are possible, although in 
many cases they are complicated as an analytic continuation from Euclidean time 
to Minkowski time is needed. For this reason most of the existing calculations have 
been performed in the quenched approximation where fermionic degrees of freedom 
are taken into account only as valence quarks and the influence of virtual quark loops 
on the gluonic background is neglected. Still these calculations prepare the ground 
for future fully dynamical calculations of transport properties of strong-interaction 
matter, for which in some cases exploratory work has already started. 

There remain many obvious, open questions in the study of the phase diagram 
of strong-interaction matter and the properties of the strongly interacting medium 
created in heavy ion collisions which can be addressed in future lattice QCD cal¬ 
culations and which will profit from future increasing computing resources as well 
as algorithmic advances. Calculations at non-zero baryon number density and the 
search for the critical point, the calculation of transport properties with fully dy¬ 
namical, light quarks as well as calculations of thermodynamic quantities using chi¬ 
ral fermion discretization schemes are among the most demanding problems which 
hopefully can be addressed in the coming years with improved techniques and re¬ 


sources. 
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